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To  my  Parents 


Preface 


This  work  on  organic  polymer-based  integrated  optic  devices  is  a  direct  extension  of 
work  done  by  Captain  Carl  A.  Kutsche  at  the  Air  Force  Institute  of  Technology  during 
1990-91.  His  investigation  of  processing  techniques  for  the  polymers  under  study  was  an 
essential  step  toward  the  design,  fabrication,  and  optical  characterization  of  working 
devices.  Captain  Kutsche's  work  concentrated  on  developing  the  processing  steps  needed 
to  fabricate  optical  waveguides  from  thin  polymer  films.  This  thesis  is  the  next  logical 
step.  The  devices  were  designed  using  the  dimensions  of  the  existing  lithographic  mask 
and  the  measured  optical  indices  of  the  polymers  under  study.  Fabrication  of  the 
waveguides  was  carried  out  using  the  techniques  developed  by  Captain  Kutsche.  Modeling 
of  the  resulting  devices  was  accomplished  using  previously  developed  approximate 
analytical  methods.  The  optical  characterization  of  the  devices  under  study  was 
accomplished  through  an  end-fired  input  coupling  technique.  The  results  of  this  thesis 
demonstrate  that  the  organic  polymers  under  investigation  have  applications  in  the  fields  of 
integrated  optics,  optical  signal  processing,  and  all-optical  computing. 

There  were,  of  course,  many  people  involved  in  this  thesis  research.  I  owe  a  debt  of 
gratitude  to  Captain  Kutsche,  then  of  WL/ELR,  for  his  advice  and  time  during  the 
fabrication  of  devices.  His  experience  in  this  area  was  invaluable  and  is  greatly 
appreciated.  Dr.  Hoa  Jiang  and  Dr.  Pete  Haaland  of  Lawrence  Associates,  WL/MLPJ, 
were  responsible  for  the  growth  of  the  thin  polymer  films  under  investigation.  I  thank 
them  for  their  help  and  cooperation.  I  would  also  like  to  thank  my  advisor,  Captain  (Dr.) 
Jeffrey  Grantham,  and  my  committee  members.  Major  (Dr.)  Paul  Ostdiek  and  Dr.  Victor 
M.  Bright.  My  conversations  with  them  and  their  editorial  comments  made  this  thesis  a  far 
better  document  than  it  would  have  been  otherwise.  Major  Richard  Bagnell's  help  in  the 
laboratory  and  contribution  of  a  format  conversion  program  were  essential  to  the  successful 


collection  of  data  and  I  owe  him  a  debt  of  gratitude.  I  would  also  like  to  thank  Chris 
O'Brien  and  Bill  Trop  for  their  help  in  wafer  preparation  and  Joe  Williams  for  his  excellent 
work  with  the  scanning  electron  microscope.  My  friend,  Eric  Baenen,  also  deserves  my 
thanks  for  making  me  laugh  and  keeping  me  company  while  at  AFIT.  His  companionship 
made  this  ordeal  much  more  tolerable  than  it  would  have  been  otherwise.  Finally,  I  would 
like  to  thank  my  wife,  Laurel;  my  parents,  Michael  and  Mary  Hamm;  our  dogs.  Bandit  and 
Noel;  and  the  cats,  Alaska  and  Doobie.  Without  their  love,  patience,  and  understanding 
this  research  would  not  have  been  accomplished. 


John  N.  Berry 
19  October  1994 


Table  of  Contents 


Preface . ii 

List  of  Figures . vi 

List  of  Tables . viii 

List  of  Symbols . ix 

Abstract . xi 

I.  Introduction . 1-1 

Background . 1-1 

Problem  Statement . 1-5 

Scope . 1-5 

Approach  and  Presentation . 1-5 

II.  Theory  of  Optical  Directional  Couplers . 2-1 

Theory  of  Dielectric  Waveguides . 2-1 

The  Ray-Optic  Approach . 2-2 

The  Electromagnetic  Development . 2-7 

Coupling  Between  Waveguides  Made  from  Linear  Media . 2-11 

Optical  Nonlinearity  in  Polymers . 2-18 

Coupling  Between  Waveguides  Made  from  Nonlinear  Media . 2-22 

III.  Design  and  Modeling  of  Waveguide  and  Coupler  Structures . 3-1 

Confinement  in  Two  Dimensions  and  Marcatili’s  Method . 3-1 

Effective  Index  Method . 3-5 

Calculation  of  Coupling  Coefficients . 3-9 

Estimation  of  the  Critical  Power  for  an  NLDC . 3-12 

Models  of  Fabricated  Coupler  Structures . 3-14 

Modeling  of  Strip-Loaded  Waveguides . 3-16 

Modeling  of  Rib  Waveguides . 3-18 

Modeling  of  Embedded  Channel  Waveguides . 3-21 

Observations  Concerning  the  Modeling  of  Waveguides . 3-21 

IV.  Fabrication  of  Waveguiding  Structures . 4-1 

Strip-Loaded  Waveguides . 4-2 

Rib  Waveguides . 4-13 

Channel  Waveguides . 4-17 


IV 


V.  Experimental  Apparatus . 5-1 

End-Fire  Coupling  Method . 5-1 

Overview  of  Experimental  Apparatus . 5-2 

Alignment  of  the  Optical  Components . 5-4 

VI.  Experimental  Procedure . 6-1 

Sample  Preparation  and  Alignment . 6-1 

VII.  Experimental  Results  and  Analysis . 7-1 

Guiding  of  Light  Through  a  Polyphenylene  Slab  Waveguide . 7-1 

Coupling  in  a  Multi-channel  Directional  Coupler . 7-4 

Lateral  Confinement  of  Light  and  Coupling  Between  Guides . 7-5 

Nonlinear  Effects . 7-10 

Waveguide  Losses . 7-12 

VIII.  Summary  and  Recommendations . 8-1 

Appendix  A:  Waveguide  Dimensions  for  Single-mode  Operation . A-l 

Appendix  B:  Slab  Waveguide  Electric  Field  Profiles . B-l 

Appendix  C:  The  Effective  Index  Method  for  Strip-Loaded  Waveguides . C-l 

Appendix  D:  The  Effective  Index  Method  for  Rib  Waveguides . D-l 

Appendix  E:  Marcatili's  Method  for  Embedded  Waveguides . E- 1 

Bibliography . Bib-1 

Vita . Vita-1 


v 


Figures 

1  - 1 .  Basic  geometry  of  an  optical  directional  coupler . 1-3 

1- 2.  The  two-dimensional  distribution  of  electric  field  within  the 

coupling  region  of  an  optical  directional  coupler . 1-4 

2  - 1 .  Radiation  of  optical  energy  into  the  cover  and  substrate  in 

a  "slab"  waveguide  structure . 2-2 

2- 2.  Path  of  a  ray  undergoing  total  internal  reflection  (TIR)  in  a 

dielectric  slab  waveguide  structure . 2-5 

2-3.  Some  possible  modes  in  a  slab  waveguide . 2-10 

2-4.  Dielectric  slab  waveguide  geometry  used  as  a  model  for  the 

development  of  the  coupled-mode  equations . 2-12 

2-5.  The  periodic  exchange  of  power  with  A/3  *  0 . 2-15 

2-6.  The  periodic  exchange  of  power  when  Aft  =  0 . 2-16 

2-7.  Dependence  of  the  power-transfer  ratio,  3,  on  the  phase  mismatch . 2-17 

2-8.  Self-defocusing  and  self-focusing  in  nonlinear  materials . 2-19 

2-9.  Uncoupled  nonlinear  waveguides . 2-20 

2-10.  The  response  of  an  NLDC  to  various  input  powers . 2-24 

2- 11.  An  NLDC  functioning  as  an  all-optical  logic  gate . 2-26 

3  - 1 .  Cross  section  of  a  guide  region  embedded  in  lower  refractive 

index  materials . 3-2 

3- 2.  A  channel  waveguide  used  for  two  dimensional  confinement  of 

optical  energy . 3-5 

3-3.  Typical  strip-loaded  and  ridge  waveguide  structures . 3-6 

3-4.  Guide  geometry  used  to  analyze  the  strip-loaded  guide  by 

Marcatili's  method . 3-8 

3-5.  Typical  coupled  waveguide  structure  and  mode  profiles . 3-10 

3-6.  The  photolithographic  mask  used  to  define  the  dimensions  of 

the  waveguide  structures . 3-15 

3-7.  Typical  mode  profiles  in  a  strip-loaded  waveguide . 3-18 


vi 


3- 8.  Electric  field  profiles  for  p  =  1  and  p  =  2  in  a  rib  waveguide . 3-20 

4-  1 .  The  photolithographic  mask  used  in  fabrication  of  waveguide  structures . 4-1 

4-2.  The  major  steps  in  the  fabrication  of  strip-loaded  waveguide  structures . 4-5 

4-3.  Side  wall  roughness  of  the  strip-loaded  waveguide  structures . 4-7 

4-4.  Two  different  techniques  for  cleaving  of  samples . 4-8 

4-5.  Resulting  sample  interfaces  for  different  cleaving  techniques . 4-9 

4-6.  Typical  SEM  photograph  used  in  the  measurement  of  film 

thicknesses  and  waveguide  features . 4-12 

4-7.  Major  steps  in  the  fabrication  of  the  rib  waveguide  structures . 4-14 

4-8.  Structures  resulting  from  rib  waveguide  fabrication  attempt . 4-15 

4-9.  The  major  steps  in  the  fabrication  of  channel  waveguides . 4-17 

4-10.  Channel  waveguides  fabricated  by  etching  silicon  dioxide  and 

deposition  of  polyphenylene . 4-19 

4- 11.  Close-up  of  the  polyphenylene  channel  waveguides . 4-20 

5  - 1 .  The  basic  geometry  of  the  end-fire  coupling  method . 5-1 

5- 2.  Test  setup  for  the  optical  characterization  of  NLDCs . 5-4 

6-  1 .  The  "sunrise"  pattern  seen  on  the  output  monitor  from  the 

diffraction  of  the  input  light . 6-2 

7-  1 .  Guiding  of  light  through  the  polyphenylene  film . 7-2 

7-2.  Measured  (a)  and  modeled  (b)  intensity  profiles  for  a 

polyphenylene  slab  waveguide . . 7-3 

7-3.  Flow  of  injected  power  in  a  multi-channel  directional  coupler . 7-4 

7  -4.  Observed  confinement  of  light  in  the  lateral  direction  and 

coupling  between  adjacent  strip-loaded  guides . 7-6 

7-5.  Observed  guiding  and  coupling  of  light  injected  into  a 

W  =  5  jim  strip-loaded  waveguide . 7-8 

7-6.  Coupling  of  light  between  adjacent  strip-loaded  waveguides . 7-10 

7- 7.  Output  intensities  for  "low"  and  "high"  peak  injected  powers . 7-11 

8-  1 .  Proposed  lithographic  mask  for  further  NLDC  research . 8-3 

vii 


Tables 


3  - 1 .  Modeling  results  for  the  strip-loaded  waveguides . 3-17 

3- 2.  Modeling  results  for  the  rib  waveguides . 3-19 

4-  1 .  Parameters  used  for  growth  of  polymer  thin-films . 4-3 

4-2.  Nominal  dimensions  of  the  strip-loaded  waveguide  layers  and  features . 4-13 


vm 


List  of  Symbols 


Symbol 


ai(z),  a2(z) 


kx>  ky,  kz 


NhN2 


nc,  ng,  ng 


Pl(z),  P2(z) 


ui(x),  u2(x) 


Definition 


functions  describing  the  amplitudes  of  the  electric  fields  in  coupled 
waveguides 


separation  distance  between  guides  in  the  coupled  slab  waveguide 
structure 


arbitrary  constant  used  in  the  normalization  of  the  electric  field  profiles 


width  of  the  guide  regions  in  the  coupled  slab  waveguide  structure 


field  parameter  that  describes  the  mode  structure  in  the  guide  region 


intensity  of  the  light  within  a  guide 


ropagation  constant  of  an  electromagnetic  wave  in  free-space 


propagation  constant  of  a  waveguide  mode  m  direction  specified  by  the 
subscript 


interaction  length  of  a  directional  coupler  that  results  in  complete  coupling 
of  the  energy  from  one  guide  to  the  other 


arameter  used  as  a  measure  of  input  power  in  Jensen's,  NLDC  model 


effective  indices  of  the  guides  in  the  coupled  slab  waveguide  structure 


index  of  the  media  around  the  guides  in  the  coupled  slab  waveguide 
structure 


refractive  index  of  the  cover  (c),  guide  (g),  or  substrate  (s)  layer  in  a  slab 
waveguide  structure 


change  in  index  due  to  nonlinear  properties  of  material  and  intensity  of 
light  propagating  within  the  material 


coefficient  that  describes  the  intensity-dependent  index  shift 


functions  describing  the  power  in  the  waveguides  of  a  directional  coupler 


injected  power  level  in  a  nonlinear  directional  coupler  that  results  in  the 
power  being  split  equally  between  the  two  output  ports 


field  parameter  that  describes  the  decay  of  the  electric  field  in  the  cover 
layer;  also,  integer  mode  number  in  x  direction  for  Marcatili' s  method 


field  parameter  that  describes  the  decay  of  the  electric  field  in  the  guide 
layer;  also,  integer  mode  number  in  y  direction  for  Marcatili' s  method 


thickness  of  the  guide  layer  in  a  slab  waveguide  structure 


functions  describing  the  spatial  distribution  of  the  modes  in  the  coupled 
slab  waveguide  structure 


width  of  waveguide  in  Marcatili's,  method 


normalized  coupler  length  in  Jensen's,  NLDC  model 


ropagation  constant  in  the  z  direction  of  the  mth  mode 


IX 


A(5 

difference  in  the  propagation  constant  between  two  waveguides  (also 
called  the  phase-mismatch  between  guides) 

£x(y) 

transverse  function  that  describes  the  profile  of  the  x-directed  electric  field 
in  the  ^-direction. 

decay  coefficients  of  the  fields  in  regions  2  and  4,  respectively,  in 
Marcatili's  method 

6m 

bounce  angle,  as  measured  from  the  interface,  of  the  mth  mode 

6c 

critical  angle  for  the  interface  between  two  media  as  measured  from  the 
interface  normal 

6cgc>  6Cgs 

critical  angle  at  the  guide/cover  (gc)  or  guide/substrate  (gs)  interfaces, 
measured  from  the  interface  normal 

K  *12,  K21 

coupling  coefficient  between  linear  waveguides 

Ao 

free-space  wavelength  of  operation 

S3,  & 

decay  coefficients  of  the  fields  in  regions  3  and  5,  respectively,  in 
Marcatili's  method 

0m 

angle  of  incidence,  as  measured  from  the  interface  normal,  for  a  ray  of  the 
mth  mode 

0gc>  0gs 

Goos-Hanchen  phase  shift  experienced  upon  reflection  from  the 
guide/cover  (gc)  or  guide/substrate  (gs)  interface. 

a  material  parameter  (chi-three)  that  describes  third-order  nonlinear 
properties 

3 

power-transfer  ratio  between  coupled  guides 

e 

symbol  used  for  the  "bar"  or  straight-through  channel  in  an  optical  coupler 

symbol  used  for  the  "cross"  or  adjacent  channel  in  an  optical  coupler 

X 


AFIT/GEO/ENP/94D-0 1 


Abstract 


This  research  was  an  investigation  into  the  suitability  of  a  recently  developed 
polymer,  polyphenylene,  as  a  material  for  integrated  optical  circuits  (IOCs).  Polymers 
show  great  promise  in  the  area  of  IOCs  because  of  material  processing  advantages, 
compatibility  with  most  existing  integrated  circuit  technology,  and  relatively  strong 
nonlinear  optical  characteristics.  This  thesis  contains  an  overview  of:  dielectric 
waveguides,  linear  and  nonlinear  directional  coupler  theory;  various  models  useful  in  the 
design  and  analysis  of  optical  waveguides;  the  fabrication  of  three  different  waveguide 
designs;  the  experimental  apparatus  and  procedure  used  to  optically  characterize  the 
waveguides;  and  the  experimental  results  of  the  characterization. 

Waveguiding  of  near  infrared  light  through  polyphenylene,  in  both  a  slab  waveguide 
and  strip-loaded  guides,  was  observed  for  the  first  time  in  this  polymer.  Coupling  of  light 
between  guides  of  a  multi-channel  directional  coupler  was  also  observed.  No  definitive 
conclusions  concerning  nonlinear  effects  are  possible  due  to  the  non-identical,  multi¬ 
channel  nature  of  the  fabricated  waveguides.  The  results  of  this  research  indicate  that 
polyphenylene  is  a  candidate  for  use  in  IOCs  and  that  the  polymer  should  be  the  topic  of 
further  research. 
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Design.  Fabrication.  Modeling,  and  Optical 
Characterization  of  Organic  Polymer  Nonlinear 

Directional  Couplers 


I.  introduction 

This  thesis  reports  on  research  that  was  conducted  on  optical  directional  couplers 
(DCs).  These  devices  are  important  integrated  optic  (10)  components  of  modem  optical 
information  processing  systems.  Recent  advances  in  organic  polymers  indicate  that  these 
materials  could  be  suited  for  use  in  10  devices  such  as  linear  and  nonlinear  directional 
couplers  (NLDCs).  Polymers  may  have  cost,  processing,  and  material  property 
advantages  over  traditional  optical  materials  used  for  waveguiding.  During  the  course  of 
this  thesis  research,  optical  DCs  were  designed,  fabricated,  modeled,  and  optically 
characterized. 

Background 

Light  has  been  used  to  communicate  information  since  prehistoric  times.  Perhaps  the 
first  instance  of  this  was  the  caveman's  use  of  campfire  light  and  shadow  puppets  to  play 
out  the  stalking  of  animals.  Other  historical  examples  abound:  signal  fires  and  lighthouses 
used  to  warn  ships  of  dangerous  coast  lines,  land-  and  ship-based  optical  telegraphs,  and 
Paul  Revere's  famous  "one  if  by  land,  two  if  by  sea"  signal  from  the  Old  North  Church 
bell  tower.  All  of  these  communication  methods  rely  on  the  modulation  of  light  to  convey 
information.  Modem  optical  communication  networks  use  this  same  concept;  information 
is  encoded  and  transmitted  using  changes  in  optical  intensity,  frequency,  or  phase 
(10:124). 

The  Information  Age  is  characterized  by  the  transmission  and  processing  of  large 
amounts  of  data.  This  flow  of  information  has  created  the  need  for  communication 
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networks  and  processing  systems  with  very  fast  data  through-put  rates.  Traditional 
communication  links  are  limited  by  the  bandwidth  available  within  the  radio  frequency 
spectrum.  Fiber  optic  data  networks  and  optical  processing  systems  have  overcome  this 
limitation  by  using  light  as  the  carrier  of  information  (1 1:2). 

The  electronic  equipment  used  to  transmit,  receive,  and  process  the  information 
carried  by  an  optical  link  is  the  bandwidth-limiting  factor  in  modem  optical  data  networks. 
With  current  electronic  equipment,  the  information  carried  by  an  optical  network  must  be 
converted  into  an  electrical  signal  before  it  can  be  processed.  After  processing,  the 
information  is  often  reconverted  into  an  optical  signal  so  that  it  can  be  retransmitted  to 
another  system.  On  a  relative  time  scale,  the  conversion  process  and  electrical  processing 
are  slow.  Great  efforts  are  being  made  to  produce  hybrid  electrical-optical  and  all-optical 
computers;  processing  systems  that  can  carry  out  operations  on  information  that  is  in  an 
optical  format.  All-optical  systems  could  be  capable  of  data  rates  on  the  order  of  1016 
bits/sec,  an  exciting  prospect  when  current  electro-optical  systems  have  achieved  data  rates 
on  the  order  of  109  bits/sec  (25:864). 

The  DC  (Fig.  1-1)  is  an  essential  element  in  an  all-optical  information  processing 
system.  A  DC  is  composed  of  two  optical  waveguides  brought  into  close  proximity.  A 
waveguiding  region,  shown  as  a  gray  line  in  the  figure,  is  simply  an  optically  suitable 
material  surrounded  by  a  material  or  materials  of  relatively  lower  refractive  index.  The 
difference  in  refractive  index  between  the  materials  leads  to  total  internal  reflection  of  the 
optical  energy  within  the  waveguiding  region. 

The  electric  field  of  the  optical  energy  traveling  though  a  waveguide  is  concentrated 
within  the  guide  material  but  some  of  the  field,  the  so-called  evanescent  tail  of  the  mode, 
extends  into  the  material  surrounding  the  waveguide.  The  adjacent  waveguide  is  placed  so 
close  to  the  "excited"  guide  (the  guide  carrying  the  light)  that  the  evanescent  tail  of  the  mode 


1-2 


overlaps  into  this  second  waveguide  (Fig.  1-2).  This  results  in  the  light  being  coupled 
from  the  excited  waveguide  into  the  adjacent  waveguide. 


Input  Port  #1 


Output  Port  #1 


©"bar"  channel 


Light  coupling  to  adjacent  guide 


Input  Port  #2 


®  "cross"  channel 


Output 
>  Light 


Fig.  1-1.  Basic  geometry  of  an  optical  directional  coupler.  Light  injected  into  the  input  of 
the  "bar"  channel  is  coupled  to  the  output  of  the  "cross"  channel  (after  25:708). 


The  amount  of  coupling  between  two  waveguides  is  a  function  of  field  confinement 
within  the  guide  material,  the  separation  distance  between  waveguides,  the  length  of  the 
interaction  region,  and  the  geometry  of  the  guides.  The  first  parameter,  the  confinement  of 
the  optical  energy,  is  a  function  of  the  difference  in  refractive  index  between  the  guide 
material  and  the  surrounding  media.  The  smaller  the  index  difference,  the  less  the  optical 
energy  is  confined  to  the  guide  material.  In  a  linear  directional  coupler,  the  index  difference 
between  the  coupler  materials  is  fixed  and  the  coupling  becomes  largely  a  function  of  the 
interaction  length  (25:266).  In  a  nonlinear  directional  coupler,  the  index  difference  between 
the  coupler  materials  is  not  fixed,  and  is  instead  a  function  of  the  intensity  of  the  light 
traveling  within  the  waveguides  of  the  coupler.  An  NLDC  works  by  changing  the 
difference  in  optical  indices  between  the  guide  material  and  the  surrounding  media.  The 
polymers  under  investigation  in  this  thesis  exhibit  this  intensity-dependent  index 
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characteristic,  also  known  as  the  optical  Kerr  effect.  The  optical  index  "seen"  by  light 
traveling  though  the  polymer  depends  on  the  intensity  of  the  beam  (4).  Thus,  as  the 
intensity  of  the  beam  is  changed,  the  index  difference  between  the  guide  material  and  the 
surrounding  media  is  changed.  This  change  in  index  modifies  the  confinement  of  the 
optical  energy  and  changes  the  amount  of  coupling  between  the  two  guides.  In  this  way, 
the  flow  of  the  light  from  the  bar  channel  to  the  cross  channel  of  a  coupler  (Fig.  1-1)  can  be 
controlled  by  the  intensity  of  the  light  itself  -  light  controlling  light  (13:1582). 


Electric  field 
strength  profile 


Regions  of  higher  index 
(waveguiding  regions) 


Region  of  lower  index 


Fig.  1-2.  The  two-dimensional  distribution  of  electric  field  within  the  coupling  region  of 
an  optical  directional  coupler.  The  two  waveguides  are  so  close  together  that  the  optical 
energy  in  the  excited  guide  overlaps  into  the  second  waveguide  region  (after  11:111). 


An  NLDC  can  function  as  an  optically  controlled  switch,  a  signal-splitting  network, 
or  as  a  logic  gate  (13:1583).  In  order  to  ease  the  transition  to  all-optical  computing  and  to 
facilitate  an  interface  with  electronic  equipment,  an  NLDC  should  be  fabricated  with 
methods  that  are  compatible  with  current  integrated  circuit  manufacturing.  Traditional 
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nonlinear  optical  materials  such  as  gallium-arsenide  (GaAs)  or  lithium -niobate  (LiNb03) 
are  not  easily  incorporated  into  standard  silicon  integrated-circuit  processing  steps  (1 1:9). 

Organic  polymer  films  have  many  advantages  over  the  materials  currently  used  in  all- 
optical  systems.  Cost  for  materials  and  processing  may  be  lower  than  they  are  for  GaAs  or 
LiNb03.  Unlike  most  nonlinear  optical  materials,  the  growth  and  processing  of  polymer 
materials  is  compatible  with  most  commonly  used  silicon  integrated  circuit  fabrication 
techniques.  Polymers  can  also  exhibit  higher  nonlinear  coefficients  than  most  commonly 
used  materials,  resulting  in  lower  optical  switching  powers  and  smaller  devices.  Finally, 
the  composition  of  organic  films  can  be  easily  altered  to  produce  the  characteristics  desired 
by  a  device  designer  (28:246;  32:778,  23). 

Problem  Statement 

This  thesis  involved  the  design,  fabrication,  modeling,  and  optical  characterization  of 
directional  couplers  fabricated  from  optically  nonlinear  organic  polymers. 

Scope 

Design  tools  such  as  the  effective  index  method  (1 1:37-42),  the  field-shadow  method 
(21:2071-2102),  and  modal  analysis  (25:248-252)  were  used  to  design  three  waveguide 
structures.  Using  the  results  of  previous  work  (18),  micron-scale  polymer  films  were 
grown  on  silicon  substrates  and  lithographically  patterned  to  produce  optical  waveguides. 
The  material  properties  and  waveguide  dimensions  were  used  to  model  the  coupling 
between  adjacent  waveguides  and  the  results  of  the  modeling  were  compared  with  the 
experimental  results.  The  coupling  characteristics  of  the  resulting  devices  were  investigated 
using  an  optical  test  setup.  Suggestions  for  further  work  in  this  area  were  considered. 

Approach  and  Presentation 

The  main  goal  of  this  research  was  to  fabricate  and  optically  characterize  directional 
couplers  that  used  polyphenylene  and  a  similar  polymer,  polythiophene,  as  waveguide 
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materials.  The  fabrication  of  devices  heavily  utilized  previous  work  (18)  and  the  modeling 
of  devices  relied  principally  on  theoretical  work  by  (11,  25, 13,  and  17). 

Chapter  II  gives  a  brief  overview  of  the  theory  of  optical  waveguides  and  directional 
couplers.  Chapter  in  details  the  modeling  of  the  coupling  structures  that  were  investigated. 
The  fabrication  of  the  polymer  films  and  devices  is  covered  in  Chapter  IV.  Chapter  V 
describes  the  setup  used  to  optically  characterize  the  waveguide  structures.  Chapter  VI 
outlines  the  procedures  used  to  align  the  system  and  characterize  the  devices.  Chapter  VII 
reports  on  the  results  of  the  experimental  investigations  and  contains  a  qualitative 
discussion  of  the  experimental  observations.  Chapter  VQI  contains  a  summary  of  this 
thesis  and  recommendations  for  further  work. 
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II.  Theory  of  Optical  Directional  Couplers 

This  chapter  covers  the  basic  theory  of  dielectric  waveguides  and  optical  directional 
couplers.  For  a  more  detailed  treatment  see,  for  example,  (1 1,  25, 17,  29,  and  13).  The 
phenomenon  of  total  internal  reflection  and  waveguiding  in  dielectric  media  is  examined. 
Next,  the  coupling  between  two  adjacent  waveguides  is  discussed.  This  discussion  is  used 
to  analyze  the  case  of  a  linear  directional  coupler.  The  concept  of  an  intensity-dependent 
index,  a  nonlinear  effect,  is  presented  and,  along  with  the  development  concerning  a  linear 
directional  coupler,  is  used  to  explain  how  a  nonlinear  directional  coupler  functions. 

Theory  of  Dielectric  Waveguides 

This  development  assumes  that  the  waveguide  and  surrounding  media  are  made  from 
lossless,  non-magnetic,  homogenous,  linear  material.  The  linear  constraint  will  be  lifted  in 
the  section  concerning  coupling  between  waveguides  made  from  nonlinear  media.  The  two 
common  approaches  to  explaining  the  phenomenon  of  waveguiding,  the  ray-optic  approach 
and  the  electromagnetic  treatment,  will  be  presented.  The  ray-optic  approach  is  a  useful 
introduction  to  the  phenomenon  of  waveguiding  because  it  provides  intuitive,  physical 
insight  into  how  light  is  confined  within  a  waveguiding  region.  In  this  approach,  light  is 
modeled  by  rays  that  bounce  off  the  media  interfaces,  obeying  simple  rules  of  reflection, 
refraction,  and  resonance.  However,  the  ray-optic  approach  does  not  provide  insight  into 
the  spatial  distribution  of  the  light's  electric  field.  This  spatial  distribution  of  the  electric 
field  is  needed  to  explain  and  model  the  coupling  that  occurs  between  adjacent  waveguides. 
The  electromagnetic  treatment  provides  this  spatial  distribution  of  the  electric  field  and  is 
essential  for  a  full  understanding  of  waveguiding  and  the  coupling  between  adjacent 
guides. 
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The  Ray-Optic  Approach 


The  ray-optic  approach  treats  light  as  a  bundle  of  rays  that  travel  through  a  medium, 
reflecting  and  refracting  at  interfaces.  In  this  way,  the  path  of  the  optical  energy  can  be 
traced  through  a  structure  like  the  one  shown  in  Fig.  2-1. 

In  Fig.  2-1,  a  simple  "slab"  waveguide  structure,  consisting  of  a  layer  of  guide 
material  with  semi-infinite  layers  of  different  materials  on  either  side,  is  shown.  The  top 
layer  is  referred  to  as  the  cover  and  has  a  refractive  index  of  nc.  The  bottom  layer  is  called 
the  substrate  and  has  a  refractive  index  of  ns,  while  the  middle  layer,  the  waveguiding 
region,  has  a  refractive  index  of  rig.  The  cover  and  substrate  are  considered  to  be  semi¬ 
infinite  in  the  y  direction  so  that  reflections  occurring  at  cover-air  and  substrate-air 
interfaces  may  be  neglected.  All  layers  have  infinite  extent  in  the  x  and  z  directions.  This 
slab  waveguide  structure  can  provide  confinement  in  one  dimension,  the  y  direction. 


x 


Fig.  2-1.  Radiation  of  optical  energy  into  the  cover  and  substrate  in  a  "slab"  waveguide 
structure.  Note  that  the  bounce  angle  associated  with  this  ray,  is  measured  from  the 
interface  while  the  critical  angles  and  angle  of  incidence,  <j>m,  are  measured  from  the  normal 
to  the  interface.  The  ray  strikes  the  cover  and  substrate  interfaces  at  an  angle  less  than  the 
critical  angles  for  the  two  interfaces.  The  cover  and  substrate  layers  have  refractive  indices 
that  are  lower  than  the  refractive  index  of  the  guide  layer  and  are  semi-infinite  in  the  y 
direction.  All  layers  are  assumed  infinite  in  the  x  and  z  directions  (after  17:35, 1 1:25). 
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For  a  ray  traveling  through  the  guide  layer  and  striking  either  the  cover  or  substrate 
layer  interface,  there  is  a  critical  angle  at  which  total  internal  reflection  (TIR) 
occurs.  The  critical  angle  is  the  smallest  angle  from  the  interface  normal  at  which  total 
internal  reflection  can  occur.  The  critical  angles  are  given  as  (25:1 1) 


QCgc=  sin  1 


(  \ 
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where  the  angles  are  measured  from  the  normals  to  the  interfaces.  Note  that  for  a  critical 
angle  to  exist,  and  thus  TIR  and  waveguiding  to  occur,  rig  must  be  greater  than  nc  and  ns. 

If  a  ray  travels  through  the  waveguide  media  and  encounters  an  interface  at  either  the 
cover  or  substrate  at  an  angle  of  incidence  less  than  the  critical  angle,  there  will  be  both 
reflection  and  refraction  of  the  ray  at  the  interface.  That  is,  some  of  the  optical  energy  will 
be  reflected  back  into  the  waveguiding  region  and  some  will  be  transmitted  through  the 
interface  and  refracted.  Fig.  2-1  shows  a  ray  striking  the  interfaces  at  an  angle  less  than  the 
critical  angles  for  the  two  interfaces.  For  a  ray  propagating  in  this  way,  some  energy  is  lost 
into  the  cover  and  substrate  during  reflections  from  the  interface.  This  type  of  ray  is  called 
a  cover  and  substrate  radiation  mode  and  will  die  out  after  propagating  a  short 
distance  down  the  guide.  The  optical  energy  is  not  guided  but  rather  is  lost  into  the  layers 
surrounding  the  guide. 

When  a  ray  is  reflected  from  an  interface,  there  will  be  a  shift  in  the  phase  of  the 
wave-front  represented  by  the  ray.  If  the  indices  of  the  cover  and  substrate  layers  are 
different  then  the  phase  shift  upon  reflection  will  be  different  for  each  interface.  The  phase 
shifts  can  be  interpreted  as  a  slight  penetration  of  the  ray  into  the  cover  and  substrate  layers. 
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A  more  detailed  discussion  of  these  shifts  in  phase,  commonly  called  the  Goos-Hanchen 
phase  shifts,  can  be  found  in  (17:27). 

Assuming  a  transverse-electric  (TE)  wave  (electric  field  perpendicular  to  the  plane  of 
incidence)  having  components  Ex>  Hy,  and  Hz,  the  phase  shifts  upon  reflection  from  the 
cover  and  substrate  interfaces  are  given  as  (11:28,  25:204) 
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where  </>gs  is  the  phase  shift  at  the  guide-substrate  interface;  (t>gc  is  the  phase  shift  at  the 
guide-cover  interface;  and  6m  is  the  bounce  angle,  as  measured  from  the  interface,  of  the 
ray  in  the  guide  medium. 

If  a  ray  strikes  an  interface  at  an  angle  of  incidence  greater  than  the  critical  angle  then 
the  ray  will  be  entirely  reflected  back  into  the  waveguide  region.  The  "zig-zag"  path  for  a 
ray  propagating  in  this  way  is  shown  in  Fig.  2-2. 

A  ray  that  undergoes  TIR  can  be  a  guided  mode  if  the  ray  also  meets  the 
transverse  resonance  condition  (16:20).  This  condition  states  that  for  a  ray  to  be 
guided  through  the  structure,  the  wave  front  must  reproduce  itself  as  it  propagates  along  the 
guide.  Mathematically  this  implies  that  the  net  phase  shift  from  a  substrate-cover-substrate 
trip  through  the  guide  medium  and  the  Goos-Hanchen  phase  shifts  from  reflections  must 
equal  a  multiple  of  2k.  The  condition  for  a  guided  mode  can  be  stated  as  (1 1 :27): 


2 k0rigtg  sin  6m  -  2 (j>gS  -  2 <j)gC  =  2mn 
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where  k0  is  the  free-space  propagation  constant  (kQ  =  oVc),  tg  is  the  thickness  of  the 
guiding  layer,  and  m  is  the  mode  number  (m  =  0,  1,  2, ...). 


cover 


Fig.  2-2.  Path  of  a  ray  undergoing  total  internal  reflection  (TIR)  in  a  dielectric  slab 
waveguide  structure.  Note  that  6m,  the  bounce  angle  associated  with  this  ray,  is  measured 
from  the  interface  while  the  critical  angles  and  angle  of  incidence,  <f>m,  are  measured  from 
the  normal  to  the  interface  (after  17:35). 


The  transverse  resonance  condition,  eqn.  (2-3),  leads  to  the  idea  of  allowed  modes 
within  the  waveguide.  Only  a  discrete  set  of  bounce  angles  will  satisfy  this  condition  and 
result  in  optical  energy  propagating  down  the  waveguide.  Rays  that  are  incident  at  angles 
greater  than  the  critical  angle  but  do  not  met  the  transverse  resonance  condition  are  lost 
through  self-interference  and  their  power  is  reflected  back  out  of  the  guide.  Rays  at  angles 
of  incidence  that  are  greater  than  the  interface  critical  angles  and  that  are  allowed  by  the 
transverse  resonance  condition  result  in  total  internal  reflection  and  the  guiding  of  energy 
along  the  z  direction.  For  a  given  set  of  nc,  rig,  ns,  tg,  and  co  (or  free-space  wavelength  of 
operation,  Ao),  the  set  of  bounce  angles  (6m,  m=  1,2, ...)  that  correspond  to  guided 

modes  can  be  calculated  using  eqns.  (2-3)  and  (2-2).  The  cut-off  condition  for  a  mode 
is  the  point  where  the  angle  of  incidence  (0m)  becomes  equal  to  the  cut-off  angle,  which  is 
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the  larger  of  the  two  critical  angles  for  the  interfaces.  A  mode  with  an  angle  of  incidence 
that  is  smaller  than  this  cut-off  angle  will  not  propagate  down  the  guide  region  and  is  said 
to  be  "cut-off." 

The  propagation  constant  in  the  z  direction  for  a  mode,  (3m,  is  given  by  (25:250) 

Pm  =  ngk0  cos(6m)  (2-4) 

Since  the  propagation  constant  of  each  mode  is  different,  a  given  mode  propagates 
down  the  guide  at  a  group  velocity  different  from  the  group  velocity  of  any  other  mode. 

All  rays  travel  through  the  guide  media  at  the  same  phase  velocity.  The  difference  in  group 
velocities  between  the  modes  is  due  to  the  fact  that  a  large  part  of  the  energy  in  higher  order 
modes  actually  travels  in  the  cover  and  substrate  layers.  One  interpretation  of  the  Goos- 
Hanchen  phase  shifts  replaces  the  shifts  in  phase  upon  reflection  with  an  additional  distance 
traveled  along  z  at  the  speed  of  light  in  the  material  across  the  interface  (either  the  cover  or 
substrate),  which  is  greater  than  the  speed  of  light  within  the  guide  region  (25:257). 
Higher-order  modes  in  a  dielectric  waveguide  undergo  greater  Goos-Hanchen  phase-shifts 
and  thus  travel  down  the  guide  faster  than  lower-order  modes  (25:255-258). 

A  waveguide  is  called  multi-mode  if  propagation  of  more  than  one  mode  is 
supported  by  the  guide.  In  a  multi-mode  guide,  a  pulse  of  light  that  excites  a  large  number 
of  modes  will  be  broadened  or  spread  in  time  as  it  travels  down  the  guide  because  the 
modes  are  not  propagating  down  the  guide  with  the  same  group  velocity.  This  pulse 
spreading  -  also  sometimes  called  modal  dispersion  -  is  due  to  the  differing  group 
velocities  of  the  modes  within  the  waveguide  (25:284-286). 

A  single-mode  waveguide  is  a  waveguide  that  allows  only  one  guided  mode  to 
propagate.  Generally,  waveguide  dimensions  must  be  small  (on  the  order  of  the 
wavelength  of  operation)  for  a  waveguide  to  be  single-mode.  In  a  single-mode  waveguide, 
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any  mode  other  than  the  zero-order  mode  (m  =  0)  results  in  an  angle  of  incidence  that  is 
less  than  the  cut-off  angle,  indicating  an  unguided  mode.  Single-mode  operation  of  a 
waveguide  has  many  advantages  in  optical  communication  systems  and  10  circuits.  The 
pulse  spreading  experienced  by  light  in  a  multi-mode  guide  can  be  a  bandwidth-limiting 
factor;  short  duration  pulses  are  "blurred"  by  the  spreading  of  the  pulses.  Single-mode 
guides,  because  of  less  pulse  spreading  and  other  effects  from  higher-order  modes,  are 
capable  of  data  transmission  rates  that  are  significantly  higher  than  the  rates  achievable  with 
multi-mode  guides  (25:286).  In  addition,  the  analysis  of  the  coupling  between  single-mode 
waveguides  is  greatly  simplified  since  the  light  propagates  down  the  guide  in  only  one 
mode.  Appendix  A  is  a  Mathcad®  (version  3.1  for  the  Macintosh,  published  by  MathSoft, 
Inc.)  file  that  uses  the  transverse  resonance  condition  to  calculate  the  waveguide 
dimensions  required  for  single-mode  operation. 

The  Electromagnetic  Development 

This  section  contains  a  brief  development  of  the  electromagnetic  equations  that  govern 
the  propagation  of  light  through  waveguiding  structures.  The  solution  of  these  equations 
results  in  a  function  that  describes  the  spatial  distribution  of  the  energy  in  the  waveguide 
structure.  As  previously  mentioned,  this  spatial  distribution  is  need  to  fully  analyze  the 
coupling  between  adjacent  waveguides.  The  section  also  discusses  how  these  equations 
lead  to  a  cut-off  condition  for  modes  within  a  waveguide.  This  development  principally 
follows  the  work  presented  in  (1 1,  34). 

The  geometry  of  the  simple  slab  waveguide,  as  shown  in  Fig.  2-1,  is  the  starting 
point  for  the  development  of  the  electromagnetic  equations  governing  the  guiding  of  light 
through  dielectric  structures.  Note  that  rig  is  assumed  to  be  greater  than  both  nc  and  ns, 
indicating  that  the  light  is  guided  by  the  structure. 

Any  traveling  electromagnetic  field  can  be  represented  by  a  combination  of  transverse 
magnetic  (TM)  and  transverse  electric  (TE)  waves.  In  this  analysis,  a  purely  TE  wave, 
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consisting  of  only  Ex,  Hy,  and  Hz  components,  is  assumed  for  simplicity.  For  such  a  TE 
wave  traveling  in  the  z  direction,  Maxwell's  wave  equation  can  be  stated  as  (11:30) 


V2Ex  = 
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df 


(2-5) 


where  tip  is  the  index  of  the  media  that  the  wave  is  propagating  through  (nc,  tig,  and  ns  for 
the  slab  waveguide  structure  of  Fig.  2-1),  and  Ex  is  the  x  component  of  the  electric  field 
propagating  in  the  z  direction.  Note  that  since  we  are  assuming  a  TE  wave,  the  Ey  and  Ez 
components  of  the  electric  field  are  zero. 

Equation  (2-5)  has  solutions  in  the  form  of  (1 1:30) 


Ex(y,z,t)  =  ex( 


(2-6) 


where  €jc(y)  is  the  transverse  function,  a  function  that  describes  the  spatial  distribution  of 
the  x  directed  electric  field  in  the  y  direction;  ©is  the  angular  frequency  of  the  radiation;  and 
(5  is  the  propagation  constant  in  the  z  direction.  The  transverse  function,  Ejy),  has  only  y 
dependence  since  we  are  assuming  that  the  waveguide  structure  is  infinite  in  the  x  and  z 
directions. 

The  transverse  function  has  the  general  form  of  (34:921) 
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where  C,  h,  q,  and  p  are  all  constants  that  describe  the  spatial  distribution  of  the 
electromagnetic  energy  within  the  waveguide  structure.  Note  that  the  field  in  the  cover  and 
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substrate  layers  decays  exponentially  because  the  optical  energy  is  concentrated  in  the  guide 
layer  for  a  guided  mode. 

Equation  (2-7)  for  the  different  regions  can  be  substituted  into  eqns.  (2-6)  and  (2-5) 
to  determine  q,  h,  and  p  (11:31,  34:921): 


h  =  (n*k02-p2f12 

(2-8b) 

q  =  (p2 -n-ck02f 2 

(2- 8a) 

P  =  (P2  -n2k02f2 

(2- 8c) 

where  k0,  previously  defined  as  oVc,  can  also  be  expressed  as  k0  =  2n/X0 ,  where  k0  is  the 
wavelength  of  the  electromagnetic  energy  in  free  space. 

The  constants  q,  h,  and  p  are  all  given  in  terms  of  p,  the  propagation  constant  in  the  z 
direction  and  the  material  constants.  Imposing  the  boundary  conditions  that  £x  and 
d£x/dy  must  be  continuous  at  the  guide-cover  and  guide-substrate  interfaces  leads  to  a 
condition  on  ft  (34:921).  In  the  resulting  transcendental  equation,  eqn.  (2-9),  both  the  left 
and  right  hand  sides  are  functions  of  /?,  k0  (and  thus  the  wavelength  of  the  electromagnetic 
energy),  and  the  material  constants  (1 1:31-32): 


q  +  p 


J 


(2-9) 


Solving  eqn.  (2-9)  results  in  a  set  of  discrete  values  of  (3  that  correspond  to  allowed 
modes  within  the  waveguide  structure.  For  each  value  of  pm,  the  values  of  pm,  hm,  and 
qm  for  each  allowed  mode  can  be  calculated  using  eqn.  (2-8).  Some  possible  modes  for  a 
slab  waveguide  structure  are  shown  in  Fig.  2-3. 
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Fig.  2-3.  Some  possible  modes  in  a  slab  waveguide.  The  profiles  depict  the  spatial 
distribution  of  the  electric  field  and  are:  (a)  a  higher-order,  non-guided  mode;  (b)  a 
substrate  radiation  mode;  (c)  the  TE2  mode;  (d)  the  TEi  mode;  (e)  the  TEo  mode.  The 
subscript  of  TE  indicates  the  value  of  m  for  that  mode  and  corresponds  to  the  number  of 
nulls  in  the  electric  field  profile  (after  11:18). 


When  a  mode  is  "cut-off,"  the  field  in  either  the  cover  or  substrate  layers  becomes 
oscillatory  instead  of  exponentially  decaying.  Cut-off  occurs  when  the  appropriate  field 
parameter,  q  or  p,  becomes  imaginary.  There  are  equations  that  can  be  used  to  find  the 
number  of  modes  that  either  a  symmetric  {nc  =  ns)  or  highly  asymmetric  (ns  »  nc) 
waveguide  will  support  (11:34-36).  Note  that  in  an  asymmetric  guide,  the  lowest  order 
mode  in  the  TE  case  can  be  cut-off.  That  is,  there  is  a  minimum  guide  layer  thickness 
required  to  support  the  lo west-order  for  a  given  set  of  materials.  This  is  in  contrast  to  the 
case  of  the  symmetric  guide  where  the  lowest  order  mode  has  no  cut-off.  However,  a 
guided  mode  will  be  poorly  confined  if  the  refractive  index  difference  between  the  guide 
and  the  cover  or  substrate  layers  is  very  small  (1 1 :34).  The  Mathcad®  document  in 
Appendix  B  can  be  used  to  model  the  profile  of  the  electric  field  in  a  slab  waveguide. 
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Coupling  Between  Waveguides  Made  from  Linear  Media 


The  basic  geometry  of  a  typical  directional  coupler  is  shown  in  Fig.  1-1.  Here,  two 
optical  waveguides  are  brought  into  close  proximity  and  allowed  to  run  parallel  to  each 
other  before  separating.  Each  of  the  waveguides  is  assumed  to  be  single-mode  and  made 
from  a  material  whose  optical  response  is  linear.  The  separation  distance  between  guides  is 
assumed  to  be  small  enough  such  that  the  evanescent  tail  of  the  field  pattern  in  one  guide 
overlaps  into  the  other  guide  (Fig.  1-2). 

The  interaction  region,  the  area  where  the  two  guides  are  close  together  and  coupling 
is  facilitated,  can  be  modeled  by  two  slab  waveguides  as  shown  in  Fig.  2-4.  The  two 
guide  regions  have  refractive  index  of  rig  and  the  surrounding  medium  has  a  refractive 
index  of  n.  For  guiding  purposes,  n  is  assumed  to  be  less  than  tig  .  The  width  of  the 
guide  regions  is  d  and  they  are  separated  by  a  distance  2a.  All  materials  are  assumed 
lossless.  The  interaction  length  needed  for  complete  transfer  of  optical  energy  from  guide  1 
to  guide  2  is  designated  as  L0  . 

Note  that  following  analysis  does  not  account  for  the  confinement  of  the  energy  in 
two  dimensions,  as  would  be  the  case  in  most  practical  10  circuits.  In  that  case,  the  spatial 
description  of  the  field  of  the  mode  would  be  a  function  of  both  the  x  and  y  coordinates. 
Even  though  the  analysis  presented  in  this  section  ignores  confinement  in  the  x  direction,  it 
is  still  illustrative  of  the  concept  of  optical  coupling  (25:264-267). 
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Fig.  2-4.  Dielectric  slab  waveguide  geometry  used  as  a  model  for  the  development  of  the 
coupled-mode  equations.  The  profiles  represent  the  distribution  of  the  electric  field.  At  zi 
the  light  is  mostly  in  guide  1,  at  Z2  the  light  is  equally  divided  between  the  two  waveguides, 
and  at  Z3  the  light  has  coupled  into  guide  2  (after  25:264). 


Coupled-mode  theory,  also  known  as  the  weak-coupling  approximation,  is  a 
common  simplified  approach  to  solving  the  problem  of  two  closely-spaced  waveguides 
(25:264).  The  theory  assumes  that  the  spatial  shape  of  the  waveguide  modes  is  not  affected 
by  the  coupling  and  that  coupling  affects  only  the  amplitudes  of  the  fields  in  the 
waveguides.  Therefore,  the  amplitudes  of  the  electric  fields  are  a  function  of  the  z 
direction,  and  are  designated  as  aj(z)  and  a2(z).  The  functions  aj(z)  and  a2(z)  describe 
how  the  optical  energy  is  distributed  between  the  two  guides  along  the  propagation 
direction,  z.  It  can  be  shown  (25:264-265)  that  the  amplitude  functions  aj(z)  and  02(1) 
obey  the  coupled-mode  equations  (25:265): 

“T'  =  ~iK 21  exp(;'A^)a2  (z)  (2-  10a) 
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and 


— =  -JK12  exp(jApz)ai(z)  (2-10b) 

dz 

where  K21  and  Kj2  are  coupling  coefficients  that  describe  how  strongly  the  modes  in  the 
two  guides  interact,  and  Aft  is  the  difference  between  the  propagation  constants  of  the  two 
guides  (also  called  the  phase  mismatch).  If  the  modes  in  the  two  guides  are  identical,  then 
A[3  =  0  and  k2i  =  Kj2  =  k.  In  general,  the  coupling  coefficients  are  (25:265) 

2 

*21  =  ^(ng  ~ n2)^j%i(y)u2(y)dy  (2-1  la) 

and 

2 

*12  =  \{nl  ~n2)^Ll-dU2('yS>Ul^dy  (2-1  lb) 

where  u](y)  and  u2(y)  are  functions  that  describe  the  transverse  profiles  of  the  modes 
within  waveguides  1  and  2,  respectively.  They  are  examples  of  the  transverse  mode  profile 
function,  eqn.  (2-7).  /)/,  are  the  propagation  constants  in  the  z  direction  of  the  modes 
within  waveguides  1  and  2,  respectively. 

The  constant,  C,  appearing  in  the  spatial  distribution  functions  uj(y)  and  U2(y)  is 
arbitrary  (1 1:32,  34:921).  In  this  development,  the  normalization  condition  for  these  two 
functions  is  taken  as  (25:253) 

EL“mW^  =  1  <2-12> 

where  um(y)  is  the  spatial  distribution  of  the  electric  field  in  the  mth  guide.  Eqn.  (2-12) 
can  be  used  to  determine  a  normalized  value  of  the  arbitrary  constant  C  in  eqn.  (2-7). 

The  phase  mismatch  per  unit  length  between  the  two  guides  is  given  by  (25:265) 
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(2-13) 


A/?  =Pi~fh. 

=  y~^Ni  -  n2) 

A0 

where  N\  =  n\  cos(#i) ;  iV2  =  n2  cos(02);  and  iV/,  N2  are  the  effective  indices  for  guides  1 
and  2,  respectively;  dj,  62  are  the  bounce  angles  of  the  modes  in  guides  1  and  2, 
respectively;  and  Xq  is  the  wavelength  of  the  electromagnetic  energy  in  free  space  (25:265). 
The  concept  of  an  effective  index  will  be  discussed  in  more  detail  in  Chapter  HI.  Again, 
for  the  special  case  of  identical  guides,  Af3  =  0  (25:266). 

From  eqns.  (2- 10a)  and  (2- 10b),  the  coupled-mode  equations,  it  can  be  seen  that  the 
variation  of  the  optical  energy  in  a  guide  depends  on  the  amount  of  optical  energy  in  the 
other  guide  and  the  coupling  coefficient.  For  the  case  of  an  initial  optical  power  launched 
into  waveguide  1  and  no  light  entering  waveguide  2  at  z  -  0,  these  equations  can  be  solved 
to  give  the  optical  power  in  the  waveguides  as  a  function  of  z,  the  direction  of  propagation 
(25:266): 


2  ( ABT  2 

Pl(z )  =  P\(0)  cos2  yz  +  —  sin2  yz 

\2Y) 


2 

P2U)  =  P1(0)^-sinV 


(2- 14a) 


(2- 14b) 


where  Pi(z),  P2U)  are  the  optical  powers  in  guides  1  and  2,  respectively;  P](0)  is  the  initial 
power  launched  into  waveguide  1  at  z  =  0;  y2  =  (A/3  /  2)2  +  k2  ;  and  K  =  -\/ki2^21  ■ 

If  the  two  guides  are  not  identical  then,  in  general,  the  phase  mismatch  will  not  be 
zero  and  the  power  will  be  exchanged  between  the  two  guides  as  shown  in  Fig.  2-5.  Note 
that  because  of  the  phase  mismatch,  the  power  exchange  is  not  complete.  That  is,  there  is 
no  interaction  length  such  that  all  of  the  power  in  guide  1  is  transferred  over  to  guide  2. 
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Fig.  2-5.  The  periodic  exchange  of  power  with  Af3  *  0.  The  period  of  the  exchange  along 
the  z  axis  is  2jt/y( after  25:266). 

If  the  guides  are  identical  then  they  have  identical  propagation  constants  and  Aft  =0. 

In  this  case  the  eqns.  (2- 14a)  and  (2- 14b)  simplify  to  (25:266) 

Pl(z)  =  Pi (0) cos2  kz  (2- 15a) 

and 

P2(z)  =  P1(0)sin2Kz  (2-15b) 

This  results  in  complete  power  transfer  between  the  two  guides.  At  z=7i/2k  the 
optical  power  in  guide  1  is  zero  and  the  optical  power  in  guide  2  is  Pi(0),  the  power  that 
was  initially  launched  into  guide  1.  This  interaction  length,  z=n/2K,  is  shown  as  L0  in 
Figs.  2-4  and  2-6.  The  periodic  nature  of  the  complete  power  exchange  between  guide  1 
and  guide  2  is  shown  in  Fig.  2-6. 
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Fig.  2-6.  The  periodic  exchange  of  power  when  Af3  =  0.  This  is  also  called  the  phase- 
matched  case  (after  25:266). 


In  the  case  of  a  directional  coupler  with  an  interaction  length  of  La,  the  power- 
transfer  ratio  can  be  defined  as  the  amount  of  power  present  in  guide  2  at  z  =  La  divided  by 
the  power  incident  in  guide  1  at  z  =  0.  This  quantity,  denoted  as  5,  is  given  by  (25:267) 


p2(  v  _ 

rn']2  .  2 
—  sine  ( 

1 

kAMl 

1/2' 

Pl(0) 

<2J 

2 

V  n  ) 

where  sinc(x)  is  defined  as  sin( tcx)/( nx). 


(2-16) 


A  plot  of  this  equation  is  shown  in  Fig.  2-7.  As  can  been  seen,  the  power- transfer 
ratio  is  a  maximum  for  Aj5  =0  (the  phase  matched  case).  Also  note  that  no  power  will  be 
transferred  to  guide  2  (3  =  0)  if  A (3L0  =  V3tt.  Thus,  the  amount  of  optical  energy  exiting 
guide  2  can  be  varied  from  Pi(0)  to  zero  by  controlling  the  phase  mismatch,  A/3 . 
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Fig.  2-7.  Dependence  of  the  power-transfer  ratio,  5,  on  the  phase  mismatch.  Note  that 
complete  coupling  occurs  only  if  the  phase-mismatch  (A(3 )  is  zero  (after  25:267). 

The  amount  of  power  transferred  from  guide  1  to  guide  2  is  a  function  of  the  guide 
materials  and  dimensions  (through  the  coupling  coefficient  k),  their  interaction  length,  and 
the  difference  in  the  propagation  constants  (phase  mismatch)  between  the  two  guides.  Any 
one  or  all  of  these  parameters  may  be  varied  to  control  the  amount  of  power  transferred 
between  the  two  guides.  Of  course,  once  a  coupler  is  fabricated  in  an  10  circuit,  the 
geometry  and  materials  of  the  coupler  are  fixed.  Note  that  a  directional  coupler  can  be  a 
"passive"  device.  Here,  passive  means  that  no  external  inputs  are  needed  to  induce 
coupling.  Optical  energy  will  transfer  from  one  guide  to  another  simply  because  the  guides 
are  in  close  proximity  and  thus  coupled. 

By  changing  the  refractive  index  of  one  or  both  of  the  guides  the  propagation  constant 
mismatch  can  be  dynamically  controlled.  In  the  case  of  lithium  niobate  or  gallium  arsenide 
structures,  for  instance,  this  mismatch  can  be  accomplished  through  the  electro-optic  effect, 
where  the  material  exhibits  a  refractive  index  that  depends  on  the  electric  field  within  the 
material  (25:268).  An  applied  voltage  can  change  the  refractive  index  in  one  or  both  of  the 
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guides  and  thus  control  the  amount  of  optical  power  that  is  transferred  from  one  guide  to 
another.  This  is  the  way  that  electro-optic  directional  couplers  route  optical  power  from 
one  output  port  to  the  other  output  port.  Thus,  an  electrically  controlled  directional  coupler 
can  function  as  either  a  switch  for  optical  energy  or  as  an  intensity  modulator.  A  directional 
coupler  made  from  nonlinear  material  functions  in  much  the  same  way.  Instead  of  an 
applied  voltage  controlling  the  coupling,  the  intensity  of  the  light  in  the  waveguides 
controls  the  amount  of  coupling  between  the  two  optical  waveguides. 

Optical  Nonlinearity  in  Polymers 

In  a  third-order  nonlinear  medium,  the  electric  field  of  an  incident  monochromatic 
light  beam  causes  a  nonlinear  polarization  to  occur  in  the  material  (24:15).  One  type  of 
third-order  nonlinearity  is  an  intensity-dependent  index  of  refraction,  also  called  the 
optical  Kerr  effect.  In  organic  polymers,  the  source  of  this  intensity-dependent  index  is 
loosely  bound  electrons  that  are  free  to  move  up  and  down  a  polymer  chain  in  the  presence 
of  an  electric  field,  such  as  the  field  of  an  incident  optical  beam  (24:35-58).  This 
polarization  of  the  atoms  or  molecules  in  the  material  can  induce  a  small  change  in  the 
optical  index  of  refraction  of  the  material.  For  more  specific  information  on  optical 
nonlinearity  in  polymers  see,  for  example,  (24). 

The  change  in  the  index  of  refraction  due  to  the  intensity  of  an  incident  beam  of  light 
is  given  by  (25:752) 


A n  =  n2I  (2-17) 

where  I  is  the  intensity  of  the  incident  light  and  «2  is  the  nonlinear  coefficient  with  units  of 
cm2/W.  This  coefficient  is  related  to  the  material  parameters  by  (25:752) 

n2  X  (2-18) 

n  e0 
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where  r\0  is  the  impedance  of  free  space,  n  is  the  refractive  index  of  the  material  with  no 
light  incident  upon  it  (the  "linear"  index),  £0  is  the  permittivity  of  free  space,  and  %<3)  is  a 
material  parameter  that  describes  the  third-order  nonlinear  properties  of  the  material 
(25:752, 24:15-19).  The  refractive  index  of  a  nonlinear  material  is  thus  a  function  of  the 
intensity  of  the  light  that  is  traveling  through  it.  This  intensity-dependent  refractive  index  is 
given  as  (25:752) 

n(I)  =  n  +  n2I  (2-19) 

Note  that  in  eqn.  (2-18),  the  sign  of  tf3)  can  be  either  negative  or  positive.  A 
negative  value  of  tf3)  indicates  that  the  refractive  index  of  the  material  decreases  as  the 
intensity  of  the  light  increases.  If  the  shape  of  the  beam  is  Gaussian  in  nature,  this  causes  a 
self-defocusing  effect  as  seen  in  Fig.  2-8a.  A  positive  value  of  %f3)  means  that  the 
refractive  index  of  the  material  increases  as  the  intensity  of  the  light  increases,  leading  to 
self-focusing  of  a  Gaussian  light  beam  as  seen  in  Fig.  2-8b. 


Fig.  2-8.  Self-defocusing  and  self-focusing  in  nonlinear  materials.  In  (a),  the  material  has 
a  negative  %f3K  The  intensity  profile  of  the  beam  creates  a  region  of  lower  refractive  index 
in  the  middle  of  the  beam  and  this  index  gradient  makes  the  beam  spread  out  and  become 
more  planar.  In  (b),  the  material  has  a  positive  The  higher  intensity  toward  the  center 
of  the  beam  creates  a  self- waveguide  region  and  a  self-focusing  effect  (after  25:755). 
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To  illustrate  the  effect  of  a  self-induced  change  in  refractive  index,  consider  two 
identical  single-mode  waveguides  made  from  a  nonlinear  material  having  a  positive  tf3) 
value.  The  guides  are  spaced  far  enough  apart  that  there  is  negligible  coupling  between 
them.  Two  separate  beams,  of  the  same  wavelength  but  different  intensities,  Ij  and  I2,  are 
launched  into  the  waveguides  at  the  same  instant  in  time.  In  this  case,  Ij  is  much  larger 
than  I 2-  At  the  entrance  to  the  guides  the  two  beams  are  in  phase  with  one  another.  This  is 
shown  in  Fig.  2-9.  It  should  be  noted  that  the  following  discussion  is  somewhat 
qualitative  in  that  the  optical  beams  propagating  through  the  waveguides  are  treated  as  plane 
waves.  A  more  complete  analysis  would  include  the  self- focusing  properties  of  the 
material,  the  effects  of  modes  in  the  waveguide  structure,  and  the  idea  of  an  effective  index 
within  the  guide  (27,  25,  24). 


Fig.  2-9.  Uncoupled  nonlinear  waveguides.  Two  beams  of  light,  with  the  same 
wavelength  and  in-phase  but  with  different  intensities  (Ij  » I2 ),  are  launched  into  identical 
nonlinear  waveguides  at  the  same  instant.  The  two  beams  become  out  of  phase  due  to  the 
intensity-dependent  refractive  index  of  the  nonlinear  material  (after  25:752-753). 


Since  the  guide  medium  is  nonlinear,  the  index  of  refraction  in  the  guide  depends  on 
the  intensity  of  the  wave  traveling  through  it.  The  guides  are  identical  but  the  intensity  of 
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the  optical  energy  traveling  through  them  is  different,  thus  each  wave  "sees"  a  different 
index  of  refraction.  These  indices  are  given  by  (25:752) 

rigi  =  n  +  n2I\  (2-20a) 


ng2=n  +  n2l2  (2-20b) 

where  ngj  is  the  refractive  index  seen  by  the  optical  energy  traveling  in  guide  1  and  ng2  is 
the  refractive  index  seen  by  the  optical  energy  traveling  in  guide  2. 

The  difference  in  index  of  refraction  between  the  two  guides  leads  to  different 
propagation  constants  for  each  of  the  waveguides.  The  mismatch  in  the  propagation 
constants  is  give  by  (25:265) 


AJ3  =ft -02 
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ih~h) 


'  (2-21) 


As  can  be  seen  in  eqn.  (2-21),  the  amount  of  mismatch  between  the  guides  is  directly 
proportional  to  the  difference  in  the  optical  intensity  traveling  through  the  guides.  Also 
note  that  if  Ij  =  I2  then  A(3  =  0  and  the  guides  are  not  phase  mismatched.  Since  the 
nonlinear  coefficient  n2  is  usually  quite  small,  the  guides  are  very  nearly  phase  matched  if 
h  =  0  and  /;  is  small  enough  not  to  induce  significant  changes  in  the  index  of  the  guide. 
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This  simplified  example  of  uncoupled  nonlinear  guides  illustrates  that  the  amount  of 
phase  mismatch  between  otherwise  identical  nonlinear  guides  can  be  controlled  simply  by 
changing  the  amount  of  optical  power  launched  into  one  of  the  guides. 

Coupling  Between  Waveguides  Made  from  Nonlinear  Media 

Previously,  the  case  of  two  nonlinear  waveguides  separated  by  a  distance  large 
enough  such  that  there  was  negligible  coupling  between  them  was  discussed.  The  phrase 
"large  enough"  implies  that  the  evanescent  tails  of  the  modes  do  not  significantly  overlap. 

In  this  section,  it  is  assumed  that  the  two  guides  have  been  brought  into  physical  proximity 
in  the  manner  illustrated  in  Fig.  2-4.  The  previously  presented  concept  of  a  passive  coupler 
and  the  phenomenon  of  an  intensity-dependent  index  of  refraction  are  brought  together  in 
this  section  to  explain  how  a  nonlinear  optical  directional  coupler  (NLDC)  functions. 
According  to  Leutheuser  (19),  this  analysis  was  first  performed  by  Jensen  (13)  and  Maier 
(20). 

It  is  assumed  that  the  waveguides  are  identical  and  single-mode.  They  are  made  from 
a  nonlinear  material  and  they  are  embedded  in  a  nonlinear  material  of  an  appropriate  index 
such  that  guiding  occurs  (13:1580).  The  analysis  of  this  geometry  leads  to  the 
development  of  four  coupling  coefficients  that  account  for  both  the  linear  and  nonlinear 
field  interactions.  Like  the  coupling  coefficients  in  the  case  of  linear  coupler,  two  of  these 
coefficients  are  a  function  of  the  guide  geometry  and  dimensions.  The  new  coefficients 
arise  from  the  nonlinear  interaction  of  the  field  with  itself  and  from  the  nonlinear  interaction 
between  the  fields  of  the  two  waveguides  (13:1581). 

The  analysis  carried  out  by  Jensen  is  involved  but  similar  in  concept  to  that  of  (25) 
for  a  linear  waveguide.  As  Jensen  points  out,  a  special  case  of  interest  for  all-optical 
computing  is  where  all  of  the  optical  power  is  initially  launched  into  just  one  waveguide 
(guide  1  in  this  example).  For  this  case,  the  power  in  guide  1  as  a  function  of  the  coupling 
length  is  (13:1582) 
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Pj(Z)  =  Pj(0){1  +  cn(2Z\mc)}  /  2 


(2-22) 


where  P/  is  the  power  in  guide  1,  Z  is  the  direction  of  propagation  normalized  by  the 
appropriate  coupling  coefficient,  Pi(0)  is  power  launched  into  guide  1,  and  cn(2Z\mc)  is  a 
Jacobi  elliptic  function  (1:573-575).  The  parameter  mc  is  a  measure  of  the  power  incident 
in  guide  1  and  is  given  by  (13:1582) 


mr  = 


[7*1(0)]" 


(2-23) 


where  Pc  is  the  critical  power. 

The  critical  power  is  defined  as  the  input  power  into  guide  1  such  that  the  power  is 
split  equally  between  the  two  guides  after  propagating  a  distance  of  Z=  k  (see  Fig.  2-10). 
This  critical  power  level  is  a  function  of  the  linear  and  nonlinear  coupling  coefficients, 
which  are  themselves  functions  of  the  material  parameters  and  guide  dimensions  (13:1582). 

The  Jacobi  elliptic  function  in  eqn.  (2-22)  is  periodic.  Therefore,  due  to  the  intensity- 
induced  phase-mismatch  at  "high"  power  levels  ( P](0 )  >  Pc ),  there  is  coupling  to  guide  2 
but  it  is  not  complete  (not  all  of  the  energy  in  guide  1  is  coupled  to  guide  2).  If  the  input 
power  is  very  small  (Pj(0) «  Pc )  then  m~  0  and  eqn.  (2-22)  can  be  approximated  as 
(13:1582): 


PX{Z)  =P1(0){l  +  cos(2Z)}/2 
=  Pi(0)cos2(Z) 


(2-24) 


Eqn.  (2-24)  is  the  form  of  the  power  exchange  for  a  linear  coupler  (13:1582,  25:266). 

The  amount  of  power  in  guide  1  as  a  function  of  distance  traveled  along  the  coupling 
region  and  mc  is  shown  in  Fig.  2-10.  The  vertical  axis  shows  the  amount  of  power  in  the 
guide  normalized  by  the  critical  power.  At  lower  power  levels,  the  guide  functions  as  a 
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linear  coupler.  That  is,  for  low  input  powers,  there  is  a  distance  ( ZJ%  ~  0.5)  where  all  of 
the  incident  power  in  guide  1  is  coupled  over  into  guide  2.  As  can  be  seen  in  Fig.  2-10, 
when  the  input  power  is  equal  to  or  greater  than  the  critical  power  the  NLDC  no  longer 
functions  in  a  manner  similar  to  a  linear  coupler.  No  matter  how  long  the  interaction 
length,  the  optical  energy  will  not  completely  couple  into  the  adjacent  guide.  At  incident 
powers  greater  than  the  critical  power,  only  a  relatively  small  percentage  of  optical  power  is 
coupled  into  guide  2  (27:101, 13:1582). 


0.0  0.5  1.0  1.5  2.0 

Z/7t 


Fig.  2-10.  The  response  of  an  NLDC  to  various  input  powers.  This  figure  illustrates  how 
the  coupler  makes  the  transition  between  acting  as  a  linear  coupler  (transferring  all  optical 
energy  over  to  guide  2  for  "low"  input  powers)  and  a  nonlinear  coupler  where  "high" 
incident  power  results  in  very  little  power  coupled  into  guide  2  for  any  interaction  length 
(after  13:1582). 


2-24 


The  work  by  Jensen  (13)  derived  the  fundamental  response  of  an  NLDC  using 
coupled-mode  theory.  The  coupling  that  occurs  between  the  waveguides  in  an  NLDC  can 
be  qualitatively  understood  in  several  different  ways.  As  seen  in  eqn.  (2-21),  the  optical 
nonlinearity  of  the  material  gives  rise  to  an  intensity-dependent  phase  mismatch  between  the 
two  guides.  This  "detuning"  (13:1582)  between  the  guides  changes  the  coupling  between 
the  guides,  as  was  seen  in  eqns.  (2-14)  for  the  case  of  a  linear  coupler.  So,  in  this  regard, 
the  NLDC  can  be  thought  of  as  a  linear  coupler  that  has  a  mechanism  for  self-modulation  of 
the  phase  mismatch.  The  intensity  of  the  light  in  the  guides  controls  the  phase  mismatch 
between  the  guides,  routing  the  optical  energy  from  one  guide  to  another  in  a  directional 
coupler  as  the  input  intensity  is  changed. 

Another  way  to  view  the  functioning  of  an  NLDC  centers  around  the  self-focusing  or 
self-defocusing  properties  of  the  nonlinear  material.  Note  that  the  coupled  mode  theory  of 
Jensen  is  based  on  the  intensity-induced  phase  mismatch  ("detuning")  between  the  guides 
and  does  not  account  for  the  self-focusing  or  defocusing  material  properties.  As  can  be 
seen  in  Fig.  2-8,  a  mode  profile  of  a  Gaussian  beam  traveling  through  a  nonlinear  material 
is  altered  by  the  intensity-dependent  refractive  index  changes  in  the  material.  For  a  self- 
defocusing  material  (a  negative  material)  the  mode  spreads  out  of  the  guide,  resulting 
in  more  overlap  with  the  adjacent  guide  and  increased  coupling.  In  a  self-focusing  material 
(a  positive  tf3)  material)  the  mode  is  more  confined  to  the  guide  through  the  self-waveguide 
effect  (Fig.  2-8).  This  increased  confinement  decreases  the  coupling  between  the  guides. 

In  this  view,  the  NLDC  acts  as  a  linear  coupler  that  has  an  intensity-dependent  coupling 
coefficient.  Thus,  the  coupling  from  one  guide  1  to  guide  2  over  a  fixed  length  depends 
upon  the  power  launched  into  guide  1. 

This  intensity-dependent  coupling  property  of  NLDCs  makes  them  suitable  for  all- 
optical  logic  devices,  such  as  AND  or  XOR  gates,  as  in  Fig.  2-11  (13:1582-1583).  The 
length  of  the  interaction  region  is  chosen  such  that  most  of  the  input  power  couples  to  guide 
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2  for  "low"  input  powers  (say  mc  =  0.1,  Z/it  ~  0.5  in  Fig.  2-10).  The  intensity  that  results 
from  the  combination  of  signals  1  and  2  when  they  are  both  "high"  is  greater  than  the 
critical  power  for  the  coupler  ( mc  >  1).  Two  simultaneously  "high"  input  signals  result  in 
detuning  of  the  guides  and  the  power  remains  in  guide  1.  However,  if  only  one  "high" 
signal  is  present,  the  signal  does  not  significantly  detune  the  guides  and  most  of  the  power 
couples  over  into  guide  2.  The  relatively  large  magnitude  and  fast  response  time  (==  10~14 
sec)  of  the  optical  nonlinearity  in  organic  polymers  makes  these  materials  strong  candidates 
for  10  devices  (13:1583,  24:2-6). 


Intensity  launched 
into  bar  channel 


Signal  #1 
1  1  1 
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Signal  #2 
1  0  1 
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Signals  #1  and  #2  are  assumed  to  be 
incoherent.  Thus,  the  result  of  their 
addition  is  simply  the  sum  of  their  intensities. 


Processed  results 
of  signals  #1  and  #2 
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Fig.  2-11.  An  NLDC  functioning  as  an  all-optical  logic  gate.  Light  pulses  of  high 
intensity,  the  result  of  simultaneously  "high"  input  signals,  exit  the  output  of  the  bar 
channel.  Pulses  of  relatively  low  intensity,  the  result  of  only  one  input  signal  being  "high," 
couple  to  the  cross  channel.  The  NLDC  thus  functions  as  all-optical  logic  element. 


Another  possible  mode  of  operation  for  the  NLDC  is  in  the  sub-critical  power  regime. 
Note  that  in  Fig.  2-10,  the  curves  for  mc  <  1  have  a  relative  shift  in  zero  points  as  the 
normalized  interaction  length,  Z,  increases.  This  is  due  to  the  nonlinear  properties  of  the 
guides  and  the  A j5Z  phase  shift  between  the  guides.  A  device  could  be  constructed  with  an 
interaction  length  such  that  complete  transfer  of  power  from  guide  1  to  guide  2  would  occur 


2-26 


at  low  input  powers.  For  this  device,  a  higher  input  power  -  but  still  less  than  the  critical 
power  -  would  result  in  a  relatively  high  percentage  of  the  input  power  exiting  guide  1. 

The  device  could  operate  as  a  logic  gate,  sensor-overload  protector,  or  signal  router  at 
powers  less  than  the  critical  power  (8). 

Recently,  work  has  been  has  been  done  to  refine  the  coupled-mode  theory  (2,  6,  22). 
The  results  of  these  investigations  indicate  that  there  are  small,  but  perhaps  in  some  cases 
significant,  differences  between  the  coupled-mode  approach  of  Jensen  and  actual  device 
performance.  Efforts  have  been  concentrated  on  accounting  for  the  self-focusing  and  de- 
focusing  effects  of  the  nonlinear  material  (6)  and  incorporating  the  change  in  mode  field 
shape  as  a  function  of  the  propagation  distance  for  the  strongly  coupled  case  (22). 
Coupled-mode  theory  assumes  that  the  interaction  between  the  fields  of  the  waveguides 
does  not  affect  the  shape  of  the  mode,  but  only  the  amplitude  of  the  electric  field.  The 
authors  of  (2)  show  that  this  is  not  always  a  good  approximation,  particularly  when  the 
spacing  between  guides  is  small  and  the  interaction  length  is  short,  as  would  be  the  case  in 
practical  10  circuits  to  be  used  for  all-optical  computing. 

Refinements  of  the  coupled-mode  theory  indicate  that  nonlinear  effect  on  the  shape  of 
the  mode  within  the  guide  region  is  very  small.  However,  outside  of  the  guide  region, 
were  the  evanescent  tail  of  the  mode  propagates,  nonlinear  effects  on  mode  shape  are 
significant  (2).  These  changes  influence  the  power  transfer  between  guides  in  a  nonlinear 
directional  coupler  because  the  coupling  is  a  function  of  the  overlap  between  the  evanescent 
tail  of  the  electric  field  and  the  adjacent  guide  region. 

The  authors  of  (2)  find  that,  in  general,  the  critical  power  calculated  when  the 
nonlinear  evanescent  tail  effects  are  accounted  for  is  lower  than  the  value  predicted  from  the 
standard  coupled-mode  approximations.  They  also  explicitly  demonstrate  that  the  small 
dimensions  required  for  practical  10  circuits  lead  to  nonlinear  effects  on  the  shapes  of  the 
evanescent  tails  of  the  waveguide  modes. 
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III.  Design  and  Modeling  of  Waveguide  and 
Coupler  Structures 

In  this  chapter,  the  design  and  modeling  of  optical  waveguides  and  directional 
couplers  are  discussed.  The  general  idea  of  confinement  in  two  directions  is  presented. 
Methods  for  estimating  the  electric  field  profiles  of  the  guided  optical  energy  in  various 
structures  of  interest  are  discussed.  The  calculation  of  the  coupling  length  between  weakly- 
coupled  waveguides  and  a  method  for  estimating  the  critical  power  of  an  NLDC  are 
presented.  Finally,  the  results  of  the  design  and  modeling  carried  out  for  this  research  are 
discussed. 

Confinement  in  Two  Dimensions  and  MarcatilP s  Method 

In  the  previous  discussions  concerning  the  theory  of  optical  waveguides,  only  the 
case  of  confinement  in  one  dimension  was  considered.  Many  practical  IO  circuits  use 
waveguides  that  confine  the  energy  in  two  dimensions.  These  optical  pathways  can  be 
used  much  like  metal  traces  on  a  standard  IC,  or  can  be  used  as  structures  in  IO  devices 
such  as  phase-modulators  or  NLDCs  (11:36,  25:239-240,  13:1582-1583).  The 
confinement  of  optical  energy  in  two  dimensions  can  be  accomplished  by  simply 
surrounding  the  guide  region  with  media  of  lesser  refractive  index  on  four  sides  instead  of 
two.  Such  a  structure  is  sometimes  called  an  embedded  guide.  The  general  form  of  this 
structure  can  been  seen  in  Fig.  3-1.  Note  that  the  view  in  Fig.  3-1  is  along  the  direction  of 
propagation,  z,  with  power  coming  out  of  the  plane  of  the  page,  as  opposed  to  a  side-view 
as  in  Fig.  2-1. 

The  media  surrounding  the  guide  region  in  Fig.  3-1  does  not  need  to  be  the  same  on 
all  four  sides  to  achieve  confinement.  The  surrounding  mediums  could  have  several 
different  indices  of  refraction,  provided  that  the  indices  of  these  surrounding  mediums  are 
all  less  than  the  index  of  the  guide  region.  The  analysis  of  this  general  structure  is 
complicated  by  the  geometry  and  boundary  conditions  that  must  be  met  at  each  interface 
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(1 1:37).  Marcatili  (21)  has  developed  an  approximate  solution  to  this  general  case  by 
assuming  that  the  modes  of  the  structure  are  well-confined  to  the  guide  region.  This 
assumption  implies  that  the  field  in  the  surrounding  regions  decays  exponentially  away 
from  the  guide  region  interfaces  and  that  the  field  in  the  comer  regions  is  negligible  (11 :38, 
21:2073).  Because  of  the  small  fields  in  the  comer  regions,  this  analysis  is  sometimes 
called  the  field-shadow  method. 


Fig.  3-1.  Cross  section  of  a  guide  region  embedded  in  lower  refractive  index  materials. 
The  view  is  along  the  direction  of  propagation  and  the  geometry  is  similar  to  that  of  the  slab 
waveguide  in  Fig.  2-1.  The  main  difference  between  the  two  guides  is  that  light  can  be 
confined  in  the  x  direction  with  this  structure  (after  1 1 :37). 


Maxwell's  equations  are  solved  by  assuming  sinusoidal  fields  in  the  guide  region  and 
exponentially  decaying  fields  in  the  surrounding  regions.  The  boundary  conditions  are 
matched  along  the  four  sides  of  the  guide  region.  Marcatili  states  that  if  the  indices  of  the 
surrounding  media  do  not  differ  by  more  than  50%  from  the  index  of  the  guide  region,  the 
modes  are  at  near  grazing  angles  and  are  essentially  transverse  electromagnetic  (TEM)  in 
nature.  This  means  that  both  the  electric  and  magnetic  fields  of  the  modes  are  transverse  to 
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the  direction  of  propagation  (9:383,  21:2074,  2096).  The  electric  field  has  two  possible 
principal  polarizations,  one  with  the  electric  field  in  the  x  direction  and  the  other  with  the 
electric  field  in  the  y  direction.  All  other  polarization  can  be  decomposed  into  combinations 
of  these  two  directions.  For  simplicity,  only  the  x  direction  polarization  is  considered  in 
this  development. 

A  TEM  wave  is  assumed  propagating  in  the  z  direction,  polarized  such  that  the  main 
transverse  field  components  are  Ex  and  Hy  (see  Fig.  3-1,  21:2075,  2078).  With  this 
polarization  and  the  geometry  of  Fig.  3-1,  the  analysis  results  in  the  following 
transcendental  equations  (21:2079) 


kxW  =  pn  -  tan  1 

(2  ) 

-tan  1 

(„2  'S 

4*^5 
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kytg=qn- tan  1  (^772) -tan  l{kyT]^j  (3-lb) 

where  the  tan'1  functions  are  to  be  taken  in  the  first  quadrant;  kx  and  ky  are  the  propagation 
constants  in  the  x  and  y  directions,  respectively;  p  and  q  are  mode  numbers  that  correspond 
to  the  number  of  peaks  in  the  field  distributions  in  the  x  and  y  directions,  respectively  (thus 
the  fundamental  mode  in  both  directions  is/?,  q  -  1);  and  where  77  are  decay  coefficients 
in  the  regions  specified  by  the  subscript  and  are  given  by 
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v  =  2,  3,  4,  5 


(3-le) 


The  transcendental  eqns.  (3-1)  are  solved  to  find  the  propagation  constants  (kx  and 
ky )  and  the  decay  coefficients  for  the  various  regions.  The  spatial  distributions  of  the  fields 
in  the  guide  and  surrounding  regions  are  described  by  these  propagation  constants  and 
decay  coefficients,  in  much  the  same  manner  as  the  analysis  of  eqns.  (2-7)  through  (2-9). 
See  Appendix  E  for  a  numerical  implementation  of  Marcatili's  analysis. 

Marcatili  states  that  this  method  is  based  on  the  assumption  that  n;  is  only  slightly 
larger  than  the  other  indices  but  that  this  assumption  is  not  very  demanding;  the  results  are 
valid  even  when  nj  is  50  percent  larger  than  112, 113,  n.4,  and  ns  (21:2074).  Comparisons  of 
this  method  with  the  full  numerical  solution  to  the  boundary  value  problem  indicate  that  the 
field-shadow  method  of  Marcatili  is  most  valid  for  aspect  ratios  (W/tg)  larger  than  two 
(21:2083).  This  is  to  be  expected  because  with  well-confined  modes,  larger  aspect  ratios 
imply  that  the  field  in  the  shadow  regions  becomes  even  less  important. 

The  channel  waveguide,  shown  in  Fig.  3-2,  is  a  special  case  of  the  general 
structure  shown  in  Fig.  3-1.  Confinement  of  the  light  in  the  x  direction  is  provided  by  the 
substrate  interfaces  in  that  direction.  Note  that  ng  =  ni,n0  =  n2,  and  ns  =  nj,  and  ny. 


Fig.  3-2.  A  channel  waveguide  used  for  two  dimensional  confinement  of  optical  energy. 
Marcatili's  method  has  been  used  to  analyze  this  structure  with  good  results  for  aspect 
ratios  (W/tg)  greater  than  two  (after  21:2083). 


3-4 


Effective  Index  Method 


It  is  possible  to  achieve  confinement  in  the  lateral  (x)  direction  without  actually 
surrounding  the  guide  region  with  materials  of  lesser  index  (1 1:41).  Confinement  is 
achieved  by  creating  a  raised  structure  on  top  of  the  guide  region,  consisting  of  either  a 
material  of  lesser  index  or  the  guide  layer  material.  Such  waveguide  structures  are  often 
called  strip-loaded  or  rib  waveguides.  Some  typical  structures  that  confine  the  light  in 
the  lateral  direction  without  embedding  the  guide  region  are  shown  in  Fig.  3-3. 

The  confinement  of  light  in  the  lateral  direction  can  be  most  easily  explained  through 
the  ray-optic  interpretation  of  waveguiding.  Consider  modes  propagating  in  regions  I  and 
II  in  Fig.  3-3a.  Note  that  the  cover  layer  in  this  case  is  not  infinite  and  thus  reflections  can 
occur  at  the  n0/nc  interface.  Because  of  these  reflections,  the  propagation  constant  of  the 
mode  traveling  in  region  I  is  less  than  the  propagation  constant  of  the  mode  in  region  II. 

The  effective  index  of  a  region  -  the  index  that  the  guide  material  seems  to  have  based  on 
the  group  velocity  of  the  mode  in  that  region  -  is  related  to  the  propagation  constant  of  the 
mode  by  (11:42) 
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where  Ne ffis  the  effective  index  of  the  region,  and  is  the  propagation  constant  in  the  z 
direction  of  the  mth  mode  in  the  region. 
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Fig.  3-3.  Typical  strip-loaded  and  ridge  waveguide  structures.  The  optical  energy  is 
traveling  in  the  z  direction  and  is  confined  to  region  II  of  the  guide  layer.  The  energy 
propagating  in  the  guide  layer  is  confined  in  the  lateral  (x)  direction  by  the  geometry  of 
either  the  cover  layer  (a)  or  the  guide  material  (b).  The  structure  in  (a)  is  a  type  of  strip- 
loaded  waveguide  while  (b)  is  referred  to  as  a  rib  waveguide  (after  1 1:42). 


Since  within  the  strip-loaded  region  (region  II)  is  greater  than  /3m  outside  the 
loaded  region  (region  I),  the  effective  index  of  the  guide  material  is  greater  under  the 
loading  strip.  This  difference  in  the  effective  indices  leads  to  confinement  in  the  lateral 
direction.  A  similar  argument  can  be  applied  to  the  rib  waveguide  structure  of  Fig.  3-3b. 

In  the  rib  waveguide,  the  effective  indices  of  the  two  regions  are  different  because  of  the 
thicker  guide  layer  in  region  II.  (11 :43). 

The  concept  of  an  effective  index  is  useful  in  analyzing  the  structures  in  Fig.  3-3. 

The  effective  index  method  was  first  introduced  by  Knox  and  Toulios  (15)  as  a  way  to 
avoid  the  explicit  matching  of  boundary  conditions  in  multi-interface  structures  such  as  Fig. 
3-3a.  The  first  step  in  the  analysis  is  finding  the  propagation  constants  in  the  z  direction  for 
the  modes  in  the  different  regions.  The  development  of  the  dispersion  relationship,  which 
relates  pm  to  the  frequency  of  operation,  geometry,  and  the  material  parameters,  is  based 
on  the  transverse  resonance  condition  and  the  Fresnel  reflection  coefficient  at  the  n0/nc 
interface.  Chuang  gives  this  dispersion  relationship  for  TE  modes  as  (5:29) 
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and  tc  is  the  thickness  of  the  cover  layer  in  the  region  (either  D  or  d  in  Fig.  3-3).  Note  that 
here  m  is  the  familiar  mode  number,  corresponding  to  the  number  of  nulls  in  the  field 
distribution.  Also  note  that  the  tan1  functions  are  to  be  taken  in  the  first  quadrant 
(15:501). 

Solving  the  transcendental  equation,  eqn.  (3-3a),  for  {5m  gives  the  propagation 
constant  in  the  z  direction  for  the  allowed  modes  in  the  region.  The  effective  index  in  the 
region  can  then  be  calculated  by  eqn.  (3-2).  Appendix  C  contains  a  numerical 
implementation  of  the  effective  index  method  applied  to  a  four-layer  waveguide  structure 
such  as  the  strip-loaded  guide  of  Fig.  3-3a.  Note  that  the  difference  in  propagation 
constants  for  the  regions  is  a  strong  function  of  the  cover  layer  thickness.  Often  D  in 
region  II  of  Fig.  3-3a  is  chosen  to  be  much  thicker  than  the  penetration  depth  of  light  into 
the  cover,  making  the  cover  layer  in  region  II  essentially  infinite.  This  condition  on  D  is 
assumed  in  the  following  equivalent  waveguide  model  (Fig.  3-4)  for  the  strip-loaded 
guides.  The  equivalent  waveguide  model  is  also  used  to  analyze  the  rib  waveguide  (Fig.  3- 
3b),  with  the  cover  region  assumed  to  infinite. 
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Fig.  3-4.  Guide  geometry  used  to  analyze  the  strip-loaded  guide  by  Marcatili's  method. 
The  equivalent  index  is  used  to  calculate  the  field  parameters  in  the  x  direction.  The  cover 
and  substrate  regions  are  assumed  to  be  semi-infinite  (after  1 1 :37) 


Equation  (3-3),  developed  for  a  four-layer  waveguide  structure,  can  also  be  used  to 
analyze  a  three-layer  waveguide  such  as  the  rib  guide  of  Fig.  3-3b.  First,  the  propagation 
constant  in  the  z  direction  is  found  for  the  two  regions  using  eqn.  (3-3).  For  this  case,  nc 
=  n0  =  1  {t]r  =  0).  To  find  the  propagation  constant  in  regions  I  and  n,  tc  in  eqn.  (3-3a)  is 
set  equal  to  the  actual  guide  layer  thickness  (either  tg  or  tg  +  D,  see  Fig.  3-3b).  The 
effective  index  method  is  used  in  Appendix  D  to  determine  the  effective  indices  for  regions 
I  and  II  of  a  rib  waveguide.  The  confinement  of  light  in  the  x  direction  can  then  be 
modeled  using  Marcatili's  method  and  an  equivalent  index  in  the  lateral  regions  around  the 
guide  region  (Fig.  3-4).  The  equivalent  index,  Neq,  used  to  model  index  of  the  lateral 
regions  for  both  the  strip-loaded  and  rib  waveguide  structures  is  given  by  (7:324) 


(3-4) 


where  Nji  is  the  effective  index  in  region  n,  and  Nj  is  the  effective  index  in  region  I. 
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Using  the  values  of  Neq,  ng,  nc,  and  ns,  the  field  distribution  parameters  for  the 
strip-loaded  and  rib  guides  can  be  calculated  using  Marcatili's  field-shadow  method 
(1 1:42).  Appendix  E  contains  a  Mathcad®  document  that  can  be  used  to  model  strip- 
loaded  and  rib  waveguide  and  the  coupling  between  waveguides  using  the  field-shadow 
method.  The  principle  advantage  in  using  the  field-shadow  method  is  that  confinement  in 
both  directions  is  accounted  for  in  determining  the  modal  field  distribution  parameters. 

Embedded  or  channel  waveguides  suffer  under  the  constraint  that  the  sidewalls  must 
be  smooth  to  reduce  losses  due  to  scattering.  In  both  strip-loaded  and  rib  waveguides, 
most  of  the  field  is  contained  in  the  guide  region,  away  from  the  sidewalls  that  define  the 
structure.  Thus,  strip-loaded  and  rib  waveguides  could  exhibit  less  scattering  losses  due  to 
side- wall  roughness  (11:44).  Excessive  roughness  and  non-uniformity  of  the  defining 
ridge  sidewalls  in  strip-loaded  or  rib  guides  could,  however,  have  a  significant  effect  on 
waveguide  losses  if  confinement  in  the  lateral  (x)  direction  is  reduced  as  a  result  of  the 
ridge  side- wall  characteristics. 

Calculation  of  Coupling  Coefficients 

The  coupling  of  light  between  slab  waveguides,  such  as  the  ones  shown  in  Fig.  2-4, 
was  examined  in  Chapter  n.  The  overlap  between  the  spatial  distributions  of  modes  in 
adjacent  waveguides  was  found  to  be  a  measure  of  the  coupling  between  the  guides.  In 
this  section,  the  special  case  of  the  coupling  between  two-dimensional,  identical,  weakly- 
coupled,  adjacent  guides  that  are  isolated  from  other  structures  is  considered  (see  Fig.  3-5). 
The  basic  process  is  the  same  as  in  Chapter  II:  the  field  distribution  parameters  for  the  two- 
dimensional  guides  are  used  to  determine  the  coupling. 

The  fabrication  process,  discussed  in  Chapter  IV,  resulted  in  groups  of  non-identical, 
perhaps  strongly-coupled,  waveguides  that  were  all  coupled  to  each  other.  This  more 
general  case  is  much  more  difficult  to  analyze  but  is,  in  principle,  similar  the  case 
considered  in  this  section;  each  guide  is  coupled  to  its  neighbor  and  the  exchange  of  energy 
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occurs  between  adjacent  guides.  Thus,  the  special  case  of  isolated  waveguides  discussed 
in  this  section  is  illustrative  of  the  phenomenon  that  would  be  observed  in  the  general  case. 
Somekh  (26)  has  analyzed  a  multi-channel  directional  coupler  consisting  of  an  infinite 
number  of  weakly-coupled  identical  guides. 


Guide  1 


Guide  2 


z 


Region  I 


Region  II 


Region  II 


Region  I 


i 


Fig.  3-5.  Typical  coupled  waveguide  structure  and  mode  profiles.  The  energy  is  traveling 
in  the  z  direction.  The  amount  of  coupling  between  the  guides  is  determined  by  the  mode 
profiles  and  their  degree  of  overlap  (after  11:111). 


For  small  input  powers,  an  NLDC  functions  as  a  linear  coupler  (13:1582).  That  is,  if 
the  intensity  of  the  light  launched  into  the  guide  is  small  enough  such  that  the  refractive 
index  of  the  guide  is  not  significantly  altered,  the  NLDC  can  be  expected  to  perform  as 
though  it  were  made  of  linear  materials.  Thus,  it  is  useful  to  calculate  the  coupling 
coefficients  of  an  NLDC  by  assuming  that  the  materials  are  linear. 

The  field  distribution  parameters  for  the  waveguide  modes  can  be  used  to  determine 
the  coupling  between  guides.  Recall  that  this  weakly-coupled  approximation  implies  that 
the  overlap  of  the  evanescent  mode-tail  in  the  guide  region  does  not  significantly  change  the 
overall  shape  of  the  mode  (11:1 16).  Marcatili  gives  the  condition  for  weakly-coupled 
guides  as  (21:2099) 
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(3-5a) 


»1 


where  s  is  the  separation  between  guides  (Fig.  3-5), 


(3-5b) 


and  kz  is  the  propagation  constant  in  the  z  direction,  which  can  be  calculated  using 
(21:2097) 


(3-5c) 


The  case  of  particular  interest  assumes  identical,  weakly-coupled,  single-mode 
waveguides  with  well  confined  modes.  Using  these  constraints,  the  coupling  coefficient 
between  guides  can  be  expressed  as  (21:2101) 


where  Kris  the  coupling  coefficient  between  adjacent  guides;  £,eq  is  the  decay  coefficient  in 
the  Neq  region  of  Figs.  3-4  and  3-5;  and  W,  s  are  the  coupler  dimensions  (see  Fig.  3-5). 
In  eqn.  (3-6a)  the  assumption  that  «2  >  ^4  has  been  made.  If  n.4  >  «2  had  been  assumed, 
rj2  in  eqn.  (3-6a)  would  have  been  replaced  by  774. 
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Equation  (3-6)  is  the  same  as  a  formula  presented  by  Marcatili  (21:2101)  except  for  a 
modification  that  accounts  for  modes  that  are  reasonably  well  confined  but  perhaps  do  not 
meet  the  strict  condition  of  eqn.  (3-5)  (26:45-46).  The  coupling  coefficient  calculated  using 
eqn.  (3-6)  is  equivalent  to  the  coupling  coefficient  presented  in  eqn.  (2-11),  provided  that 
the  proper  normalization  of  the  field  profiles  is  accomplished  using  eqn.  (2-12).  Good 
agreement  between  the  two  methods  for  calculating  k  were  observed  during  modeling  of 
the  waveguide  structures  (see  Appendix  E,  p.  E-8). 

For  identical  guides  and  modes,  the  length  required  for  complete  coupling  of  the  tight 
from  one  guide  to  the  other,  La,  is  related  to  the  coupling  coefficient  by  (1 1:1 15) 


Recall  that  the  power  will  be  transferred  periodically  between  the  two  guides,  as  indicated 
by  eqn.  (2-15)  and  Fig.  2-6. 

Estimation  of  the  Critical  Power  for  an  NLDC 

The  critical  power  for  an  NLDC  is  a  function  of  the  nonlinearity  of  the  guide  material 
and  the  mode  profile  (13:1581).  It  is  assumed  here  that  the  interaction  length  of  the  coupler 
is  L0  =  n /2  k,  where  zr  is  the  linear  coupling  coefficient  between  the  guides.  By 
approximating  the  mode  in  the  guide  as  a  plane  wave,  an  estimate  of  the  critical  power  for 
an  NLDC  can  be  obtained  as  (13:1583) 


Pc  =  kA- 


2n^n2esu 


(3-8) 


where  Kris  the  linear  coupling  coefficient  between  adjacent  guides,  A  is  the  cross-sectional 
area  of  a  guide,  ng  is  the  linear  index  of  the  guide,  c  is  the  speed  of  tight  in  vacuum,  and 
n2esu  is  the  nonlinear  coefficient  in  electrostatic  units  (esu). 
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From  the  relationship  (in  cgs-Gaussian  units),  n  =  rig  +  H2esu  I -El2,  where  E  is  the 
electric  field  in  statvolt/cm,  it  can  been  seen  that  fi2esu  must  have  units  of  cm2/statvolt2. 
These  units  are  equivalent  to  cm-sec2/gram.  That  is,  [n2esu]  =  cm2/statvolt2  =  cm- 
sec2/gram  (3:307). 

The  value  of  ri2esu  can  be  found  from  the  intensity-dependent  nonlinear  coefficient 
n2mks  (in  mks  units)  by  (3:308) 


n2 esu  ~ 


n2mksng 

8.4 


(3-9) 


where  ri2mks  must  be  expressed  in  cm^/kW  and  ri2esu  carries  units  of  cm-sec^/gram. 

The  critical  power  estimated  with  eqn.  (3-8)  is  based  on  the  assumption  that  the 
waveguide  mode  can  be  approximated  by  a  plane  wave.  Thus,  the  estimate  is  only  valid 
for  the  fundamental  mode  and  should  be  regarded  as  a  means  of  determining  the 
approximate  power  level  at  which  nonlinear  optical  effects  could  be  observed  in  an  NLDC. 

Note  that  the  assumption  was  made  that  the  interaction  length  between  the  two  guides 
was  L0  =  tz/2k.  Thus  for  very  weakly-coupled  guides,  where  Kris  very  small,  the 
interaction  length  between  guides  can  be  very  long.  Since  L0  is  large,  the  accumulated 
phase  shift,  Aj3L0,  can  be  substantial  even  for  small  values  of  intensity-induced  A (5.  Thus, 
only  small  input  powers  are  needed  to  achieve  switching  in  a  coupler  with  a  very  large 
interaction  length.  It  should  be  noted,  however,  that  coupler  structures  with  interaction 
lengths  longer  than  a  few  millimeters  are  not  practical  in  IOCs.  The  critical  powers  for 
these  structures  are  essentially  meaningless  because  the  coupling  between  the  guides  is  too 
weak  to  achieve  switching  within  reasonable  interaction  lengths. 


Models  of  Fabricated  Coupler  Structures 

Three  different  waveguide  structures,  strip-loaded,  rib,  and  embedded  channel,  were 
designed  and  modeled  (see  Figs.  3-2  and  3-3).  The  goal  was  to  design  single-mode 
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waveguides  and  to  model  the  coupling  that  occurs  between  adjacent  guides.  All  modeling 
was  conducted  at  the  nominal  wavelength  of  0.92  pm.  Known  parameters  such  as  index 
and  nonlinearity  for  the  materials  of  interest  were  used  in  the  design  process.  The  two 
polymers  used  in  the  fabrication  process,  polythiophene  and  polyphenylene,  have  indices 
of  1.71  and  1.79,  respectively,  at  0.85  pm  (14).  The  nonlinear  coefficient  of 
polyphenylene,  n2mks ,  is  2.03  x  lO'1^  cm2/kW  (n2esu  =  4.33  x  10'H  cm-sec^/gram)  at 
0.85  pm  (14).  The  absorption  resonance  of  polyphenylene  is  in  the  UV  so  the  material 
parameters  are  approximately  constant  across  the  near-infrared  spectrum  (14).  Note  that 
polyphenylene  is  a  self-defocusing  material  so  the  sign  of  the  nonlinear  coefficient  is 
negative.  The  loss  coefficient  of  polyphenylene  at  0.85  pm  has  been  measured  to  be  2.3  ± 
0.6  dB/cm  (14).  For  simplicity,  all  materials  were  modeled  as  lossless. 

An  additional  design  constraint  was  the  use  of  a  pre-existing  photolithographic  mask. 
The  dimensions  of  this  mask,  used  in  previous  work  (18),  defined  the  feature  sizes  of  the 
waveguides.  This  mask  consisted  simply  of  three  groups  of  parallel  lines,  each  group 
having  the  same  basic  features  but  with  the  separation  between  the  features  differing  from 
group  to  group.  A  diagram  of  the  2"  chromium  mask  used  to  define  the  waveguiding 
structures  is  shown  in  Fig.  3-6. 

In  the  case  of  the  strip-loaded  structures,  the  use  of  negative  photoresist  resulted  in 
waveguides  with  nominal  dimensions  of  W  =  3,  5,  and  10  pm  and  separation  distances  of 
j=l,2,  and  3  pm.  Using  positive  photoresist  in  the  fabrication  of  the  rib  waveguide 
structures  resulted  in  nominal  dimensions  of  W  -  1,  2,  and  3  pm  and  5  =  3,  5,  and  10  pm. 
For  the  channel  guide  structures,  waveguide  regions  with  nominal  dimensions  of  W  =  3,  5, 
and  10  pm  and  5=1,2,  and  3  pm  were  produced  by  using  positive  photoresist.  More 
information  concerning  the  fabrication  of  waveguide  structures  is  contained  in  Chapter  IV. 
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Fig.  3-6.  The  photolithographic  mask  used  to  define  the  dimensions  of  the  waveguide 
structures.  Note  that  two  different  sets  of  waveguide  structures  were  produced  with  this 
mask,  differing  by  the  type  of  photoresist  (positive  or  negative)  used  in  the  fabrication 
process.  The  groups  of  guides  were  separated  by  approximately  2  mm. 


The  first  step  in  the  design  process  was  the  determination  of  the  allowable  dimensions 
for  single-mode  operation.  In  the  y  direction  the  structures  were  approximated  by 
equivalent  slab  waveguides  and  the  Mathcad®  document  in  Appendix  A  was  used  to  find 
the  minimum  and  maximum  allowable  dimensions  for  single-mode  operation.  The  field- 
shadow  method  was  used  to  determine  how  many  modes  a  structure  would  support  in  the  x 
direction  by  varying  the  parameter  p,  the  mode  number  in  the  x  direction.  The  highest 
value  of  p  for  which  a  solution  to  eqn.  (3- la),  bounded  by  the  condition  that  kx  <  k0, 
exists  was  the  number  of  modes  that  the  structure  can  support  in  the  x  direction. 

The  profiles  of  the  electric  field  were  modeled  using  the  field  distribution  parameters 
calculated  by  the  field-shadow  method  of  Marcatili.  For  simplicity,  only  the  x  direction 
polarization  of  the  TEM  modes  is  considered  (see  Fig.  3-1).  The  Mathcad®  documents 
found  in  Appendices  C,  D,  and  E  were  used  to  compute  and  plot  these  profiles.  The  field 
distribution  parameters  were  also  used  to  compute  the  coupling  coefficient  and  coupling 
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length  for  adjacent  guides.  The  critical  power,  subject  to  the  constraints  on  eqn.  (3-8),  was 
calculated  where  applicable. 

The  numerical  aperture  of  the  guides  was  calculated  using  the  indices  of  the  guide 
material,  cover,  and  substrate  layers,  as  appropriate.  Since  the  guides  were  asymmetric  in 
the  x  direction,  the  index  used  in  the  calculation  was  the  larger  of  either  the  cover  or 
substrate.  This  is  because  the  larger  index,  resulting  in  a  larger  critical  angle  at  that 
interface,  is  the  limiting  factor  in  determining  if  a  ray  will  be  guided.  The  numerical 
aperture,  NA,  was  calculated  using  (25:251) 

WA=("«-"«)1/2  (3'10) 
where  nc,s  is  the  index  of  either  the  cover  or  substrate,  whichever  is  larger. 

Modeling  of  Strip-Loaded  Waveguides 

A  waveguide  structure  like  the  one  shown  in  Fig.  3-3a  was  designed  and  modeled 
using  both  the  effective-index  and  field-shadow  methods.  In  this  structure,  polythiophene 
was  used  as  the  cover  layer  ( nc  =1.71),  polyphenylene  was  used  as  the  guide  material  (tig 
=  1.79),  and  silicon  dioxide  was  used  as  the  substrate  layer  ( ns  =  1.5).  The  minimum  and 
maximum  allowable  values  for  tg  were  calculated  by  assuming  that  the  cover  and  substrate 
layers  were  semi-infinite.  It  was  found  that  the  minimum  guide  thickness  needed  to  avoid 
cutting-off  the  fundamental  mode  was  0.28  pm  while  the  maximum  allowable  thickness  for 
single-mode  operation  was  1.15  pm. 

Based  on  post-fabrication  inspection  of  SEM  photographs,  nominal  values  of  the 
dimensions  D,  d,  and  tg  were  found  to  be  1.2  pm,  0.6  pm,  and  0.67  pm,  respectively. 
Using  Appendix  E,  it  was  found  that  the  W  =  10  pm  waveguides  were  capable  of 
supporting  both  the  fundamental  (p  =  1)  and  first  odd  mode  (p  =  2)  in  the  x  direction.  The 
modeling  results  for  the  various  waveguide  widths  and  separations  are  shown  in  Table  3-1. 
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An  example  of  the  field  distribution  in  W=  10  pm  guides  is  shown  in  Fig.  3-7.  Only 
coupling  between  the  fundamental  modes  is  shown  in  this  figure.  Likewise,  only  coupling 
between  identical  modes  is  shown  in  Table  3-1.  In  general,  there  would  also  be  coupling 
between  the  fundamental  mode  in  one  guide  and  the  first  odd  mode,  if  supported,  in  an 
adjacent  guide. 


Table  3-1.  Modeling  results  for  the  strip-loaded  waveguides.  The  model  parameters  were: 
X0  =  0.92  |im,  D  =  1.2  pm,  d  -  0.6  pm,  tg  =  0.67  pm,  NA  =  0.53,  ri2esu  =  4.33  x  10'1  * 
cm-sec^/gram,  Neff  outside  guide  region  was  1.743,  Neff  within  the  guide  region  was 
1.744,  and  Neq  used  in  Marcatili's  method  was  1.789. 
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tValue  is  for  the  first  odd  mode  (p  =  2)  in  the  x  direction.  Note  that  no  critical  powers  are 
given  for  the  first  odd  mode  because  the  approximation  used  to  develop  eqn.  (3-8)  is  only 
valid  for  the  fundamental  mode.  Critical  powers  are  not  given  for  couplers  with  L0  longer 
than  a  few  millimeters  because  they  would  not  be  practical  devices. 
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Fig.  3-7.  Typical  mode  profiles  in  a  strip-loaded  waveguide.  The  thickness  of  the  guide 
material  is  0.67  |im,  the  width  of  the  guides  is  10  pm,  and  the  separation  between  the 
guides  is  3  pm.  The  zero  point  on  the  horizontal  axis  is  the  center  of  the  guide  in  that 
direction.  In  (a),  the  spatial  distribution  of  the  electric  field  in  the  y  direction  is  shown. 
Note  that  the  mode  is  asymmetric  about  the  center  of  the  guide  because  the  indices  of  the 
cover  and  substrate  layers  are  significantly  different.  The  distribution  of  the  electric  fields 
in  the  x  direction  and  their  overlap  in  a  directional  coupler  structure  are  shown  in  (b). 


Modeling  of  Rib  Waveguides 


Rib  waveguides  like  the  one  shown  in  Fig.  3-3b  were  designed  and  modeled  using 
the  effective-index  method  and  the  field  shadow  method.  The  material  parameters  were 
taken  as  n0  =  1  (air),  rig  =  1.71  (poly thiophene),  and  ns  =  1.5  (silicon  dioxide).  For  this 
structure,  it  was  found  that  the  minimum  guide  thickness  needed  to  avoid  cutting-off  the 
fundamental  mode  was  0.17  pm  while  the  maximum  allowable  thickness  for  single- mode 
operation  was  0.73  pm.  The  guide  thickness  was  chosen  to  be  0.5  pm  and,  based  on 
modeling  by  the  effective  index  method,  D  for  design  purposes  was  chosen  to  be  0.1  pm. 

The  number  of  supported  modes  in  the  x  direction  for  a  given  waveguide  structure 
was  found  using  the  field-shadow  method.  The  W  -  1  pm  guides  were  single-mode  in  the 
x  direction  while  both  the  W  =  2  and  W  =  3  pm  were  found  to  support  the  first  odd  mode. 
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Post-fabrication  inspection  of  the  rib  waveguides  indicated  that  the  polythiophene 
guide  layer  was  not  grown  thick  enough  to  produce  practical  waveguides.  The  fabrication 
process  and  inspection  results  are  discussed  further  in  Chapter  IV.  The  modeling  results 
presented  in  Table  3-2  were  calculated  using  the  design  values  for  the  guides. 


Table  3-2.  Modeling  results  for  the  rib  waveguides.  The  model  parameters  were:  XQ  = 
0.92  (im,  tg  =  0.5  pm,  D  =  0.1  }im,  NA  =  0.82,  n2esu  =  4.33  x  10"  H  cm-sec^/gram, 
Neff  outside  guide  region  was  1.612,  Neff  within  the  guide  region  was  1.632,  and  Nen  = 
1.691. 
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tValue  is  for  the  first  odd  mode  (p  =  2)  in  the  x  direction.  Note  that  no  critical  powers  are 
given  for  the  first  odd  mode  because  the  approximation  used  to  develop  eqn.  (3-8)  is  only 
valid  for  the  fundamental  mode.  Critical  powers  are  not  given  for  couplers  with  L0  longer 
than  a  few  millimeters  because  they  would  not  be  practical  devices. 


A  coupling  length,  L0,  longer  than  several  millimeters  indicates  that  the  guides  are 
essentially  uncoupled  (see  Table  3-2).  The  critical  powers  in  Table  3-2  were  calculated 
using  the  value  of  n2esu  f°r  polyphenylene  and  are  given  for  the  fundamental  mode  of 
reasonably  well-coupled  guides.  Note  that  for  a  given  geometry  the  fundamental  mode  (p 
=  1)  in  the  x  direction  can  be  uncoupled  from  adjacent  waveguides  while  the  first  odd  mode 
(p  =  2)  can  be  reasonably  well  coupled.  For  instance,  with  IT  =  2  pm  and  s  =  10  pm  the 
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fundamental  mode  has  a  coupling  length  of  L0  ~  1  x  10^  mm  while  the  first  odd  mode  has 
a  much  more  practical  value  of  L0  =  3.6  mm.  This  great  disparity  in  coupling  lengths  is 
because  the  fundamental  mode  is  very  well  confined  for  the  mentioned  geometry  but  the 
first  odd  mode  is  nearly  cut-off  and  thus  poorly  confined.  Such  a  configuration  could  be 
used  to  "strip"  the  first  odd  mode  from  a  waveguide,  leaving  only  the  fundamental  to 
propagate.  Plots  of  the  electric  field  profiles  for  the  p  =  1  and  p-  2  modes  illustrate  the 
difference  in  coupling  for  the  two  modes  (see  Fig.  3-8). 


x  (pm) 


(a) 


(b) 

Fig.  3-8.  Electric  field  profiles  for  p  =  1  and  p  =  2  in  a  rib  waveguide.  The  width  of  a 
guide  is  W  =  2  pm,  the  separation  between  guides  is  s  =  10  pm.  Guide  1  is  centered  on  x 
=  0.  Note  that  the  p  =  1  mode  (a)  is  essentially  uncoupled  while  there  is  a  large  degree  of 
overlap  (and  thus  coupling)  between  modes  for  the  p  =  2  mode  (b). 
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Modeling  of  Embedded  Channel  Waveguides 


Channel  waveguides  like  the  one  shown  in  Fig.  3-2  were  designed  and  modeled 
using  n0  =  1  (air),  tig  =  1.79  (polyphenylene),  and  ns  =  1.5  (silicon  dioxide).  The 

structure  was  analyzed  by  using  Marcatili' s  field-shadow  method.  The  thickness  of  the 
guiding  layer,  based  on  post-fabrication  inspection,  was  1.8  pm. 

Modeling  indicates  that  the  guides  can  support  multiple  modes  in  both  the  x  and  y 
directions.  For  instance,  with  W  =  10  (lm  and  tg  =  1.7  jlm,  the  guide  can  support  4  TE 
modes  in  the  y  direction  and  22  TE  modes  in  the  x  direction.  All  but  the  very  highest 
modes  are  essentially  uncoupled,  even  for  s  =  1  pm.  Both  of  these  results,  multi-mode 
guides  and  lack  of  coupling,  are  due  to  the  relatively  large  difference  in  index  between  the 
substrate  and  guide  materials.  Using  Appendices  A  and  E,  the  maximum  dimensions  for 
single-mode  operation  are  tg  ~  0.6  |im  and  W  ~  0.47  pm.  These  dimensions  imply  that 
coupling  light  into  single-mode  guides  of  this  geometry  would  be  very  difficult. 

Observations  Concerning  the  Modeling  of  Waveguides 

The  results  of  the  modeling  for  strip-loaded,  rib,  and  channel  waveguides  indicate 
that  the  strip-loaded  geometry  allows  the  largest  dimensions  for  single-mode  operation. 

The  presence  of  a  cover  layer  with  a  refractive  index  close  to  the  guide  material  refractive 
index  permits  the  single-mode  guide  thickness  to  be  larger  than  it  is  in  the  more 
asymmetric  rib  or  channel  guides.  The  strip-loaded  geometry  with  dJD  ~  0.5  weakly 
confines  the  energy  in  the  lateral  direction.  The  net  result  of  this  weak  confinement  is  that 
the  guide  can  be  relatively  large  in  the  lateral  direction  and  still  remain  single-mode.  Also, 
coupling  between  guides  is  enhanced  by  the  weak  confinement.  The  "large"  single-mode 
dimensions  of  strip-loaded  guides  make  them  good  candidates  for  experimental  work.  The 
coupling  of  light  into  a  guide  becomes  more  feasible  as  the  dimensions  of  the  guide  are 
increased. 


3-21 


IV.  Fabrication  of  Waveguiding  Structures 

This  chapter  details  the  fabrication  of  the  three  different  waveguide  structures 
designed  and  modeled  during  this  research.  Some  basic  understanding  of  integrated  circuit 
fabrication  techniques  is  assumed  in  this  chapter;  see,  for  example,  (12,  30).  Each 
waveguide  structure  was  fabricated  using  similar  wafer  preparation,  polymer  growth, 
feature  etching,  and  cleaving  techniques.  In  addition,  the  features  of  all  three  waveguide 
structures  were  defined  using  the  same  photolithographic  mask.  The  film  growth  and 
feature  etching  techniques  were  developed  during  previous  work  on  this  topic  (18). 

The  photolithographic  mask  used  for  waveguide  fabrication  was  developed  during 
previous  work  by  Kutsche  (18:32).  The  mask  was  reused  in  this  work  to  expedite  the 
fabrication  process.  As  mentioned  in  Chapter  in,  this  mask  consisted  simply  of  three 
groups  of  parallel  lines,  each  group  having  the  same  basic  features  but  with  the  separation 
between  the  features  differing  from  group  to  group.  A  diagram  of  the  2"  chromium  mask 
used  to  define  the  waveguiding  structures  is  shown  in  Fig.  4-1. 
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Fig.  4-1.  The  photolithographic  mask  used  in  fabrication  of  waveguide  structures.  Note 
that  two  different  sets  of  waveguide  structures  were  produced  with  this  mask,  differing  by 
the  type  of  photoresist  (positive  or  negative)  used  in  the  fabrication  process.  The  groups 
were  separated  by  approximately  2  mm. 
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During  the  photoresist  exposure  step,  the  mask  was  exposed  three  times  across  the 
surface  of  the  wafer,  resulting  in  three  copies  of  the  same  pattern.  This  was  possible 
because  the  mask  was  a  dark-field  mask;  it  blocked  the  light  in  the  outlying  regions,  leaving 
them  unexposed.  The  pattern  was  repeated  across  the  surface  of  the  wafer  to  maximize  the 
number  of  waveguides  produced. 

All  scanning  electron  microscope  (SEM)  photographs  were  taken  at  relatively  low 
energy  (typically  one  or  two  keV)  to  avoid  damaging  the  polymer  films.  In  addition,  the 
samples  were  coated  with  approximately  30  A  of  tungsten  to  prevent  localized  charging  of 
the  dielectric  materials  and  enhance  resolution  of  the  photographs.  All  SEM  pictures  were 
taken  with  a  Hitachi  S-900  Low  Voltage  High  Resolution  (LVHR)  scanning  electron 
microscope,  an  instrument  well  suited  for  producing  high-resolution  images  under  low- 
energy  conditions. 

Strip-Loaded  Waveguides 

The  strip-loaded  waveguides  consisted  of  two  polymer  films  deposited  on  a  silicon 
wafer  that  had  been  thermally  oxidized  to  form  the  substrate  layer  (see  Fig.  3-3a).  These 
waveguides  were  fabricated  using  negative  photoresist  resulting  in  dimensions  of  W  =  3,  5 
and  10  (im  and  5  =  1,2,  and  3  [im.  This  section  details  the  steps  used  to  fabricate  the  strip- 
loaded  guides  used  in  this  research. 

The  first  step  in  the  fabrication  process  was  the  preparation  of  the  silicon  wafers.  A 
group  of  3"  wafers  were  cleaned  by  immersing  them  in  a  3:2  solution  of  sulfuric  acid 
(H2SO4)  and  hydrogen  peroxide  (H2O2)  for  25  minutes.  They  were  then  rinsed  for  20 
minutes  in  a  flowing  deionized  (DI)  water  bath  and  blown  dry  with  nitrogen  gas.  The 
wafers  were  then  thermally  oxidized  in  a  dry  oxygen  atmosphere  for  120  hours  at  1 100°C. 
The  wafers  were  slowly  inserted  (about  1  in/min)  into  the  oxidation  furnace  to  avoid 
thermal  stress-induced  warpage  of  the  wafers.  The  oxidation  process  resulted  in 
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approximately  2.2  |im  of  high  quality  "dry"  silicon  dioxide  (SiC>2)  being  grown  on  the 
surface.  This  group  of  oxidized  wafers  were  used  as  substrates  for  the  strip-loaded  guides. 

Immediately  prior  to  film  growth,  the  wafer  to  be  used  was  re-cleaned  by  immersing 
it  in  the  3:2  H2S04:H202  cleaning  solution  for  25  minutes,  rinsing  it  in  DI  water  for  25 
minutes,  blowing  it  dry  with  nitrogen,  then  throughly  drying  it  in  a  165°C  convection  oven 
for  25  minutes.  To  avoid  contamination  of  the  polymer  film  the  wafer  was  kept  in  a 
silicon  wafer  carrier  until  it  was  inserted  in  the  polymer  growth  chamber. 

The  polymer  films  were  deposited  through  plasma-enhanced  chemical  vapor 
deposition  (CVD)  (18:8-16,  31).  Polyphenylene  was  grown  as  the  guide  layer  with 
polythiophene  serving  as  the  cover  layer.  After  the  wafer  to  be  grown  on  was  inserted,  the 
deposition  chamber  was  allowed  to  pump-down  for  three  days  prior  to  film  growth  to 
remove  any  traces  of  water  vapor  or  other  out-gassed  contaminates.  The  conditions  for 
growth  of  the  polymer  films  were  as  shown  in  Table  4-1.  The  substrate  temperature  was 
not  controlled  during  the  polymer  growth  process.  Polythiophene  was  grown  as  the  cover 
layer  by  changing  the  monomer  gas,  replacing  the  benzene  flow  with  dichlorothiophene 
(18:10-15). 


Table  4-1.  Parameters  used  for  growth  of  polymer  thin-films. 


Growth  Parameter 

Value 

Total  chamber  pressure 

0.5  torr 

Flow  rate  of  argon 

100  cc/min 

Flow  rate  of  monomer  gast 

0.6  cc/min 

Forward  RF  power 

100  W 

Reflected  RF  power 

84  W 

tBenzene  for  polyphenylene,  dichlorothiophene  for  polythiophene. 
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During  film  growth  the  wafer  was  rotated  at  approximately  4  rpm  for  more  uniform 
deposition  across  the  surface  of  the  wafer.  The  films  were  grown  without  any  sort  of 
diffuser  on  the  gas  outlet  nozzle.  Polyphenylene  was  grown  for  28  hours  and 
polythiophene  was  grown  for  55  hours.  Post-growth  inspection  of  this  deposition 
indicated  growth  rates  of  0.024  (iim/hr  for  polyphenylene  and  0.022  (im/hr  for 
polythiophene.  The  resulting  polymer  coated  wafer  was  stored  at  room  temperature  in  a 
dry  box  (<  20%  relative  humidity)  to  avoid  foreign  object  contamination  and  delamination 
from  water  vapor  (18:36). 

The  major  steps  in  the  fabrication  of  the  strip-loaded  guides  are  diagrammed  in  Fig. 
4-2.  The  first  step  in  the  patterning  of  the  polythiophene  cover  layer  was  the  plasma 
deposition  of  silicon  nitride  (Si3N4)  as  a  passivating  layer  (18:27-29).  The  growth 
conditions  for  the  nitride  layer  were  as  follows:  a  gas  mixture  of  5%  SiF4  (silane)  and 
95%  N2  at  a  flow  rate  of  168  sccm/min;  a  separate  flow  of  N2  at  1250  sccm/min;  a  flow  of 
NH3  (ammonia)  at  10  sccm/min;  forward  RF  power  of  21  W;  substrate  temperature  of 
120°C;  and  a  total  pressure  of  approximately  840  millitorr.  Growth  was  monitored  by  also 
growing  nitride  on  an  otherwise  bare  wafer  and  performing  post-growth  ellipsometry.  A 
total  of  0.12  (im  of  Si3N4  was  deposited  with  an  observed  growth  rate  of  0.01  pm/min. 

After  nitride  growth,  photoresist  was  spun  onto  the  wafer  to  prepare  for 
photolithographic  masking.  Waycoat  28  negative  photoresist  was  spun  on  to  a  thickness  of 
approximately  1  pm  and  soft-baked  at  100°C  for  5  minutes.  The  photoresist  was  then 
exposed  for  8  sec  using  the  photolithographic  mask  described  in  Fig.  4-1.  The  pattern  was 
repeated  three  times  across  the  wafer  surface  as  previously  noted.  Note  that  there  was  no 
intimate  contact  between  the  mask  and  the  wafer  due  to  the  fact  that  the  mask  was  designed 
to  be  used  with  2"  wafers.  This  accounts  for  the  rough  side- wall  features  of  the  resulting 
waveguides  (Fig.  4-3).  The  photoresist  was  developed  for  25  seconds,  rinsed  with 
isopropyl  alcohol,  and  the  wafer  dried  for  3  minutes  at  100°C  with  a  contact  hot  plate. 
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Fig.  4-2.  The  major  steps  in  the  fabrication  of  strip-loaded  waveguide  structures.  In  (a) 
the  substrate  is  shown  coated  with  polymer  films,  the  silicon  nitride  passivating  layer,  and 
negative  photoresist.  The  photoresist  layer  is  exposed  and  developed  (b)  and  then  the 
silicon  nitride  is  etched  (c).  The  polythiophene  layer  is  then  partially  removed  with  an 
oxygen  etch  (d).  Note  that  this  step  also  etches  away  the  remaining  photoresist. 
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The  wafer  was  then  placed  in  the  reactive-ion  etch  (RIE)  chamber  for  partial  removal 
of  the  cover  layer.  The  nitride  layer  was  removed  (where  exposed)  with  a  mixture  of  O2 
and  Freon- 14  known  commercially  as  DE-100.  The  nitride  was  etched  for  1.2  minutes  at  a 
rate  of  0.1  |j,m/min,  completely  removing  the  nitride  not  protected  by  the  photoresist  and 
leaving  the  polythiophene  exposed.  Exposed  polythiophene  was  etched  using  oxygen  with 
a  flow  of  50  seem  at  a  total  pressure  of  50  millitorr.  Note  that  the  photoresist  remaining  on 
top  of  the  protected  nitride  was  also  etched  during  the  removal  of  the  polythiophene.  The 
RTF,  was  capacitively-coupled  with  a  DC  bias  of  300  VDC  and  forward  power  of  98  W. 
Polythiophene  was  etched  for  a  total  of  7.5  minutes  with  an  observed  etch  rate  of  0.8 
(im/min.  This  agrees  well  with  previous  work  (18:47).  The  sample  was  removed  several 
times  during  the  polythiophene  etch  for  profilometry  as  a  check  on  the  etch  rate.  Note  that 
the  layer  of  silicon  nitride  remaining  on  the  top  of  the  waveguide  structures  (Fig.  4-2d) 
could  have  been  removed  with  another  DE-100  etch  step  but  was  left  to  provide  some 
measure  of  protection  against  damage  by  subsequent  handling.  This  layer  is  thin  enough 
and  far  enough  away  from  the  guiding  layer  to  be  ignored  in  the  modeling  of  the 
waveguides. 
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Fig.  4-3.  Side  wall  roughness  of  the  strip-loaded  waveguide  structures.  This  top  view  of 
the  sample  shows  the  very  rough  side  walls  in  the  cover  layer  of  the  waveguides.  Both  5 
pm  and  3  pm  wide  structures  are  shown.  The  separation  between  waveguides  in  this 
group  is  nominally  2  pm. 


The  wafer  was  cleaved  in  order  to  take  cross-section  pictures  of  the  resulting 
waveguide  structures.  The  samples  were  cleaved  both  in  compression  and  in  tension. 
"Compression"  and  "tension"  refer  to  the  type  of  force  applied  to  the  polymer  film.  In 
compression  cleaving,  the  back  of  the  wafer  is  scored  and  broken  by  applying  force  to  the 
silicon  side  of  the  wafer,  with  the  polymers  facing  down.  Two  slightly  different  methods 
of  tension  cleaving  were  tried.  In  one  case,  the  wafer  was  scored  across  the  back  (on  bare 
silicon  side)  and  in  the  other  a  small  tick  mark  was  made  on  the  polymer  side  of  the  wafer, 
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near  the  edge.  In  tension  cleaving  methods,  the  sample  was  broken  with  the  polymer  films 
facing  up.  The  different  cleaving  techniques  are  shown  in  Fig.  4-4.  A  lens  tissue  was 
used  between  the  glass  slide  and  the  polymer  film  to  lessen  the  damage  to  the  film. 
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Fig.  4-4.  Two  different  techniques  for  cleaving  of  samples.  In  (a),  the  polymer  film  is 
compressed  while  it  is  broken  while  in  (b)  the  film  is  cleaved  in  tension.  Note  that  two 
different  scoring  techniques  were  tried  with  tension  cleaving. 


Typical  cleaved  surfaces  produced  by  the  different  cleaving  techniques  are  shown  in 
Fig.  4-5.  The  sample  in  Fig.  4-5a  was  broken  by  the  compression  technique  shown  in 
Fig.  4-4a.  In  Fig.  4- 5b,  the  sample  was  broken  in  tension  after  being  scribed  across  the 
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sample  on  the  bare  silicon  side.  Fig.  4-5c  shows  the  results  of  cleaving  the  sample  in 
tension  after  scoring  a  small  "tick"  mark  near  the  edge  on  the  polymer  side  of  the  wafer. 


(a) 

Fig.  4-5.  Resulting  sample  interfaces  for  different  cleaving  techniques.  In  (a),  the  sample 
was  cleaved  in  compression  after  being  fully  scored  on  bare  silicon  surface.  A  sample 
cleaved  in  tension  after  being  fully  scored  on  the  bare  silicon  surface  is  shown  in  (b).  In 
(c),  a  sample  cleaved  in  tension  after  being  scored  with  a  small  "tick"  mark  on  the  polymer 
side  of  the  wafer,  near  the  edge,  is  shown.  The  interface  quality  in  (a)  and  (b)  are 
approximately  the  same  while  the  interface  in  (c)  shows  much  more  film  damage  and 
interface  roughness. 


The  goal  of  the  cleaving  process  is  to  produce  a  smooth,  clean,  perpendicular  edge  so 
that  light  my  be  efficiently  coupled  into  the  waveguide.  The  most  successful  technique 
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seemed  to  be  cleaving  in  compression  after  scoring  the  back  surface  of  the  wafer.  Breaking 
the  wafer  with  the  polymer  film  in  tension  after  fully  scoring  the  bare  silicon  side  of  the 
wafer  also  produced  acceptable  results.  Observation  of  samples  cleaved  with  these  two 
techniques  did  not  show  marked  differences  in  cleave  quality.  Both  resulted  in  relatively 
smooth  interfaces  for  the  coupling  of  light.  Cleaving  the  wafer  in  tension  after  scribing 
only  a  small  "tick"  mark  near  the  edge  of  the  wafer  resulted  in  a  large  degree  of  film 
damage  and  interface  roughness  (see  Fig.  4-5c).  This  type  of  interface  damage  was  typical 
across  the  sample  cleaved  in  tension  by  use  of  the  "tick"  mark  method.  The  differences  in 
the  cleave  quality  can  be  seen  clearly  by  comparing  Figs.  4-5a  and  4- 5b  with  Fig.  4- 5c. 


(4-5b) 
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(4- 5c) 


The  physical  size  of  the  different  layers  and  features  of  the  strip-loaded  waveguide 
structure  were  measured  from  SEM  photographs.  A  typical  picture  used  for  making  these 
measurements  is  shown  in  Fig.  4-6  and  nominal  dimensions  for  the  guide  features  are 
shown  in  Table  4-2.  No  wide-spread  delamination  of  the  polyphenylene  film  was 
observed  during  SEM  inspection  of  the  samples.  The  small  areas  of  film  damage  that  were 
observed  seemed  to  be  caused  by  the  cleaving  and  handling  processes,  not  water  vapor  or 
processing  chemicals  as  has  been  observed  in  poly  thiophene  (18). 
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The  uniformity  of  the  guiding  layer  thickness  across  the  wafer  was  determined  by 
measuring  the  thickness  of  this  layer  on  two  different  samples,  one  taken  from  near  the 
center  of  the  wafer  and  the  other  taken  ~  1.3  cm  radially  away  from  the  center.  The 
variation  of  the  film  thickness  was  observed  to  be  approximately  0.07  (am  between  the  two 
samples,  with  the  film  thickness  decreasing  away  from  the  center  of  the  wafer.  Because 
the  roughness  of  the  side  walls  precluded  accurate  measurement  of  the  separation  between 
guides  and  guide  width,  the  nominal  design  values  were  taken  for  these  dimensions. 
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Fig.  4-6.  Typical  SEM  photograph  used  in  the  measurement  of  film  thicknesses  and 
waveguide  features. 


Table  4-2.  Nominal  dimensions  of  the  strip-loaded  waveguide  layers  and  features. 


Parameter 

Nominal  Value  fuml 

0.67 

d 

0.6 

D 

1.2 

ts 

2.2 

Rib  Waveguides 

The  attempted  fabrication  of  rib  waveguides  was  not  successful  due  to  insufficient 
polymer  film  growth.  The  wafer  used  for  this  process  was  taken  from  a  stock  of  polymer 
coated  wafers  kept  by  the  growers,  WL/MLPJ.  The  wafer  was  coated  with  polythiophene 
and  had  a  film  thickness  of  ~  0.2  4m.  This  wafer  was  processed  in  much  the  same  way  as 
the  wafer  used  in  the  fabrication  of  strip-loaded  guides.  The  two  major  exceptions  were 
that  the  wafer  was  2"  in  diameter  -  thus  better  suited  for  use  with  the  2"  photolithographic 
mask  -  and  the  processing  was  done  with  positive  photoresist.  The  use  of  positive 
photoresist  resulted  in  nominal  guide  widths  of  W=  1,  2,  and  3  (im  and  separation 
distances  of  s  =  3,  5  and  10  fim.  The  major  steps  in  the  processing  of  this  wafer  are 
shown  in  Fig.  4-7. 
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Fig.  4-7.  Major  steps  in  the  fabrication  of  the  rib  waveguide  structures.  In  (a),  the  silicon 
substrate  is  shown  coated  with  polythiophene,  silicon  nitride,  and  photoresist.  The 
positive  photoresist  is  exposed  and  developed  (b)  and  the  wafer  is  etched  to  remove  the 
nitride  layer  where  exposed  (c).  An  oxygen  etch  is  then  carried  out  (d)  to  partially  remove 
the  polythiophene  layer. 


Post-processing  inspection  of  the  wafer  surface  using  the  SEM  showed  that  the  film 
thickness  was  too  thin  to  fabricate  practical  waveguides.  In  addition,  delamination  of  the 
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polythiophene  film  was  observed  in  several  locations  on  the  wafer.  Typical  features  of  the 
post-processing  structures  are  shown  in  Fig.  4-8. 


(a) 

Fig.  4-8.  Structures  resulting  from  rib  waveguide  fabrication  attempt.  This  is  a  top  view 
of  the  rib  waveguide  structures.  The  W  =  2  (im  waveguide  group  is  shown  in  (a).  A 
close-up  of  the  rib  waveguide  structure  and  film  delamination  is  shown  in  (b).  Note  how 
much  smoother  and  uniform  the  sidewalls  of  the  structure  are  as  compared  with  Fig.  4-3. 
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(4-8b) 

Fig.  4-8  shows  that  contact  between  the  mask  and  wafer  during  the  lithographic 
process  results  in  waveguide  side  walls  that  are  much  smoother  than  those  obtained  during 
the  fabrication  of  the  strip-loaded  waveguides  (see  Fig.  4-3,  note  the  difference  in  scale 
between  Figs.  4-8  and  4-3).  This  demonstrates  that  the  RIE  process  is  capable  of 
producing  smooth,  uniform  micron-scale  waveguide  structures  in  organic  polymer 
materials.  The  nature  of  the  small  "droplets"  or  spots  on  the  film  is  unknown.  They  could 
be  foreign  object  contamination  or  perhaps  spalling  of  the  material  due  to  processing  steps. 
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Channel  Waveguides 


Channel  waveguides,  similar  to  the  one  shown  in  Fig.  3-2,  were  fabricated  using 
silicon  dioxide  as  the  substrate  and  polyphenylene  for  the  guide  material.  Channels  were 
etched  into  an  oxidized  2"  wafer  using  an  RIE  process.  Polyphenylene  was  then  grown 
over  the  wafer,  filling  in  the  channels  and  forming  waveguides.  The  major  steps  in  the 
fabrication  of  the  channel  waveguides  are  shown  in  Fig.  4-9.  Note  that  the  polyphenylene 
was  also  deposited  on  the  barriers,  making  the  resulting  surface  non-planar. 
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Fig.  4-9.  The  major  steps  in  the  fabrication  of  channel  waveguides.  In  (a),  the  oxidized 
wafer  is  shown  coated  with  photoresist.  The  positive  photoresist  was  exposed  using  the 
lithographic  mask  and  then  developed  (b).  Channels  were  etched  into  the  silicon  dioxide 
(c)  and  then  polyphenylene  was  deposited  to  form  waveguides  (d).  The  waveguides  are 
the  regions  between  the  silicon  dioxide  sidewalls. 
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The  wafer  was  processed  using  positive  photoresist,  resulting  in  etched  channels  with 
W=  3,  5  and  10  pm  and  separating  sidewalls  of  silicon  dioxide  with  s  =  1,  2  or  3  pm 
wide.  Some  undercutting  during  the  etching  of  the  silicon  dioxide  changed  these 
dimensions,  generally  decreasing  the  separation  between  waveguides  and  increasing  the 
widths  of  the  guides.  The  photoresist  remaining  on  the  wafer  after  the  etching  of  the 
silicon  dioxide  was  removed  in  a  plasma  photoresist-ashing  chamber.  The  wafer  was  then 
cleaned  by  immersing  it  in  a  3:2  H2S04:H2C>2  solution  for  25,  rinsing  it  in  DI  water  for  25, 
blowing  it  dry  with  nitrogen,  and  then  throughly  drying  it  in  a  165  °C  convection  oven 
before  film  growth  began. 

Polyphenylene  was  grown  on  the  wafer  surface  with  a  plasma  deposition  process, 
resulting  in  structures  as  shown  in  Fig.  4- 10.  The  growth  conditions  for  this  film  were  the 
same  as  the  polyphenylene  grown  for  the  strip-loaded  waveguides  (shown  in  Table  4-1). 
Post-growth  inspection  of  this  film  indicated  a  deposition  rate  of  0.02  pm/hr  and  that  a  total 
of  1.8  pm  of  film  was  deposited. 

The  relative  roughness  of  the  deposited  film  and  the  poor  cleave  interface  can  been 
seen  in  Fig.  4-11.  The  sample  was  cleaved  in  tension  by  scoring  fully  across  the  bare 
silicon  (see  Fig.  4-4b).  The  roughness  of  the  silicon  dioxide  can  be  attributed  to  the  RIE 
process.  The  polyphenylene  film  deposited  during  the  fabrication  of  the  channel 
waveguides  appears  to  be  physically  different  from  the  film  grown  for  the  strip-loaded 
waveguides.  The  film  growth  for  the  strip-loaded  guides  resulted  in  what  appears  to  be  a 
dense,  smooth,  uniform  polyphenylene  layer.  As  can  be  clearly  seen  in  Fig.  4-11,  the 
polymer  film  in  the  channel  waveguides  appears  to  be  much  less  uniform  and  smooth. 

This  could  be  due  to  the  roughness  of  the  silicon  dioxide  substrate  layer  in  the  channel 
guides;  the  first  layer  of  polymer  follows  the  contours  of  the  rough  silicon  dioxide  and  that 
roughness  is  propagated  out  as  the  film  grows.  Another  possible  cause  is  contamination, 
or  changes  in  composition  due  to  time,  of  the  benzene  used  as  the  monomer  gas.  The 
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history  of  the  distilled  benzene  used  for  film  growth  in  these  two  cases  is  unknown; 
perhaps  the  benzene  used  for  growth  of  the  channel  waveguides  film  was  not  freshly 
distilled  or  contaminated. 


181314  2,8K 


Fig.  4-10.  Channel  waveguides  fabricated  by  etching  silicon  dioxide  and  deposition  of 
polyphenylene.  The  waveguides  are  the  flat  portions  between  the  raised  ridges  of  polymer. 
The  ridges  are  an  artifact  of  the  polymer  deposition  process. 
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Fig.  4-11.  Close-up  of  the  polyphenylene  channel  waveguides.  Note  the  relative 
roughness  of  the  top  of  the  waveguide  region  and  the  poor  cleave  interface. 
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V.  Experimental  Apparatus 

This  chapter  contains  the  details  concerning  the  experimental  apparatus  used  to 
optically  characterize  the  waveguides  fabricated  during  this  research.  Brief  overviews  of 
the  techniques  used  for  coupling  light  into  a  waveguide  and  aligning  the  system  are  given. 
Some  criteria  for  component  selection  are  discussed. 

End-Fire  Coupling  Method 

The  end-fire  coupling  method  was  chosen  as  the  technique  for  launching  light  into  a 
waveguide  because  of  its  flexibility,  relative  simplicity,  and  because  it  is  well  suited  for 
injecting  light  into  an  individual  waveguide  (1 1:90-92).  The  basic  geometry  of  the  end-fire 
coupling  method  can  be  seen  in  Fig.  5-1. 


Laser 


Fig.  5-1.  The  basic  geometry  of  the  end-fire  coupling  method.  On  the  input  side  light 
source  is  focused  down  onto  the  guide  layer.  Light  is  coupled  into  the  guide  and  exits  the 
other  side  where  it  is  gathered  and  imaged  by  the  output  lens  (after  1 1:91). 


For  efficient  coupling,  the  spatial  profile  of  the  input  laser  beam  should  match  the 
shape  of  the  waveguide  mode  that  is  to  be  excited.  In  addition,  the  spot  size  of  the  focused 
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beam  should  closely  match  the  waveguide  dimensions.  Finally,  the  numerical  aperture 
(NA)  of  the  input  lens  and  the  waveguide  should  be  matched,  with  the  NA  of  the  input  lens 
(NAf)  equal  to  or  less  than  the  NA  of  the  waveguide  (NA g).  If  NAi  >  NAg,  some  of  the 
light  entering  the  guide  will  strike  the  interfaces  at  less  than  the  critical  angle  and  be  lost.  If 
all  the  above  conditions  are  met  and  the  system  is  in  perfect  alignment,  the  input  coupling 
efficiency  could  approach  100%  (11:92).  In  practice,  however,  lower  efficiencies  are 
expected  due  to  mismatches  in  beam  profile,  spot  size,  NA,  and  the  difficulties  of  aligning 
a  system  with  sub-micron  tolerances.  Calculations  assuming  a  Gaussian  input  beam  with  a 
waist  of  2  (im  indicated  that  input  coupling  efficiency  could  be  as  high  as  50%  for  the  strip- 
loaded  guides  under  investigation.  After  losses  in  the  optical  beam  train  and  any  sample 
misalignment  are  taken  into  account,  30%  might  be  a  more  reasonable  estimate  for  input 
coupling  efficiency. 

The  spot  size  for  a  diffraction-limited  input  lens  can  be  calculated  from  (25:131) 

ds=  1.22^2-  (5-1) 

s  NA 

where  ds  is  the  diameter  of  the  spot  and  NA  is  the  numerical  aperture  of  the  input  lens. 
Approximately  84%  of  the  beam  power  is  contained  within  the  area  of  this  spot. 

Overview  of  Experimental  Apparatus 

The  equipment  used  to  optically  characterize  the  fabricated  devices  was  assembled 
and  aligned  on  an  optical  bench.  Fig.  5-2  is  a  top-view  diagram  of  that  equipment.  A 
wavelength  tunable  argon-pumped  Ti:Sapphire  laser  (ATS),  with  an  output  wavelength 
ranging  from  =  0.78  |im  to  ~  0.96  (im,  was  used  as  the  excitation  source.  Optical 
characterization  was  typically  performed  at  an  output  wavelength  of  0.92  jim.  In  mode- 
locked  configuration,  the  laser  was  capable  of  2  ps  wide  pulses  with  a  time  between  pulses 
of  13.3  ns.  Peak  powers  in  the  GW  range  were  obtainable  in  this  mode.  Neutral  density 


5-2 


filters  were  used  where  appropriate  to  regulate  the  intensity  of  the  beam  (NDW,  NDF).  An 
optical  isolator  (01)  was  used  to  prevent  feedback  into  the  laser  from  the  optical  setup,  and 
a  half-wave  plate  (POL)  was  used  to  reorientate  the  polarization  of  the  beam  to  the  x 
direction.  An  acousto-optic  modulator  (AOM)  was  used  to  divert  some  of  the  light  to 
another  experiment  on  the  bench.  The  zero-order  beam  was  used  in  this  setup.  Various 
irises  (I)  and  an  alignment  plate  (AP)  were  used  to  align  the  optics  to  the  beam. 

The  sample  under  test  (SUT)  was  mounted  on  a  five-axis  stage  such  that  the  laser 
beam  excited  the  guide  as  shown  in  Fig.  5-1.  The  input  microscope  objective  (MO)  had  a 
magnification  power  of  40x  and  a  numerical  aperture  of  NA{  =  0.65.  Using  eqn.  (5-1)  and 
an  operating  wavelength  of  X0  =  0.92  pm,  the  spot  size  of  this  objective  was  »  2  pm.  The 
output  objective  was  a  20x  unit  with  an  NA  of  NA0  =  0.4.  The  strip  loaded  guides  had  an 
NA  of  NAg  =  0.53.  Thus,  there  were  some  injection  losses  incurred  because  NA[  >  NAg. 
An  input  objective  with  an  NA  larger  than  NA g  was  chosen  to  decrease  the  input  spot  size. 
Note  also  that  since  NAg  >  NA0,  not  all  of  the  light  coming  from  the  waveguides  was 
gathered  by  the  output  objective.  An  output  objective  with  a  smaller  NA  was  chosen  to 
increase  the  field  of  view  on  the  output  side  of  the  SUT. 

The  images  of  the  SUT  interfaces  were  observed  with  cameras  sensitive  to  the  near 
infrared  spectrum  and  monitors.  A  personal  computer  was  used  to  capture  and  digitize  the 
output  of  the  cameras.  A  red  helium-neon  laser  (HNL)  was  aligned  to  the  optical  system 
and  used  for  the  initial  alignment  and  imaging  of  the  SUT. 
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Fig.  5-2.  Test  setup  for  the  optical  characterization  of  NLDCs.  ATS  is  an  argon-pumped 
tunable  Ti:Sapphire  laser;  P  is  a  periscope  arrangement  of  mirrors;  01  is  an  optical  isolator; 
POL  is  a  half-wave  plate  used  to  rotate  the  polarization  of  the  beam;  NDW  is  a  neutral 
density  filter  wheel;  AOM  is  an  acousto-optic  modulator;  M  is  a  flat  mirror;  I  is  an  iris  used 
in  the  alignment  process;  PBS  is  a  pellical  beam  splitter;  AP  is  an  alignment  plate  for 
aligning  the  optics  to  the  beam;  PM  is  the  power  meter  head;  MO  is  a  microscope  objective 
mounted  on  a  five-axis  stage;  SUT  is  the  device-under-test  mounted  on  a  five-axis  stage; 
BB  is  a  beam  block;  HNL  is  a  helium-neon  laser;  NDF1  and  NDF2  are  sets  of  neutral 
density  filters;  Cl  and  C2  are  cameras  for  imaging  the  input  and  output  side  of  the  SUT, 
respectively;  TV  is  a  monitor  for  observing  a  camera  output;  and  PC  is  a  personal  computer 
used  for  frame-grabbing  of  the  camera  output. 


Alignment  of  the  Optical  Components 

The  alignment  of  the  system  was  complicated  by  the  number  of  components  and  the 
sub-micron  dimensions  of  the  guiding  layer.  The  input  microscope  objective  was  aligned 
with  the  retro-reflection  technique  using  the  AP  and  an  iris  in  the  SUT  position.  The  AP 
consisted  of  a  polished  aluminum  plate  with  several  small  holes.  A  hole  in  the  AP  was 
aligned  with  the  center  of  the  input  beam  and  the  retro-reflection  of  the  input  objective  was 
observed  on  the  AP.  The  input  objective  was  aligned  such  that  the  retro-reflections  of  both 
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the  front  and  rear  surfaces  of  the  objective  coincided  with  the  hole  in  the  AP.  The  iris  in  the 
SUT  position  was  used  to  check  this  alignment  by  changing  the  z  position  of  the  objective 
and  observing  the  spot  on  the  iris.  The  output  objective  was  aligned  to  the  beam  and  the 
input  objective  using  an  iris. 

With  the  pellical  beam  splitter  (PBS)  in  place  on  the  output  side,  the  HNL  was 
aligned  such  that  its  beam  retraced  the  ATS  beam.  The  HNL  beam  was  used  to  initially 
image  the  output  side  of  the  sample  during  installation  of  the  SUT.  After  this  initial 
alignment,  the  PBS  could  be  replaced  with  a  mirror  (M)  to  increase  the  amount  of  the  light 
reaching  the  output  camera. 
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VI.  Experimental  Procedure 

This  chapter  details  the  procedures  used  to  characterize  a  sample  under  test  (SUT). 
The  preparation  of  samples  for  testing  and  the  alignment  of  the  sample  to  the  input  beam  are 
discussed.  Capture  of  images  and  simple  analysis  of  those  images  using  a  personal 
computer  is  explained.  Only  the  strip-loaded  waveguides  were  optically  characterized  due 
to  the  poor  fabrication  results  of  the  other  waveguide  designs. 

Sample  Preparation  and  Alignment 

Samples  for  optical  characterization  were  cleaved  from  the  wafer  using  the 
compression  technique  illustrated  in  Fig.  4-4.  The  length  of  the  samples,  the  z  dimension, 
ranged  from  2  to  5  mm.  Sample  preparation  met  with  varying  levels  of  success.  Finding 
a  sample  with  good  interface  cleave  quality  and  waveguides  unspoiled  by  handling  damage 
or  fabrication  flaws  along  the  length  of  the  guides  was,  at  times,  a  tedious  process.  Poor 
guide  quality  may  have  been  due  to  a  number  of  factors:  excessive  losses  due  to  the 
roughness  of  the  side  walls  (see  Fig.  4-3);  the  fact  that  each  group  of  guides  was 
concentrated  in  an  approximately  100  |im  wide  region  and  thus  susceptible  to  "mass 
destruction"  by  handling  damage;  or  the  practice  of  placing  the  polymer  side  down  during 
the  cleaving  process,  thereby  causing  film  and  waveguide  damage. 

After  cleaving,  a  sample  was  mounted  onto  a  fixture  with  a  small  amount  of  rubber 
cement.  The  rubber  cement  was  applied  to  the  fixture  surface  and  the  bare  silicon  side  of 
the  sample  was  pressed  into  the  cement,  taking  care  not  to  break  the  wafer,  contaminate  the 
interface  surfaces  with  cement,  or  scrape  across  a  region  of  the  wafer  containing  guides. 
The  fixture  was  mounted  to  a  five-axis  stage  for  alignment  of  the  SUT  and  input  beam. 

The  SUT  was  placed  into  the  beam  train  with  the  input  beam  blocked  and  the  helium 
neon  laser  (HNL)  on.  All  alignments  were  conducted  with  the  input  beam  reduced  in 
intensity  by  the  neutral  density  filter  wheels.  In  addition,  neutral  density  filters  were  used 
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as  appropriate  to  reduce  the  intensity  of  the  images  at  the  cameras.  The  output  side  of  the 
SUT  was  imaged  using  the  output  objective  and  the  HNL  beam.  This  step  insured  that  an 
image  of  the  SUT  output  edge  appeared  in  the  plane  of  the  output  camera.  The  HNL  could 
then  be  turned  off.  The  input  beam  was  focused  on  the  input  side  of  the  wafer.  The  image 
from  this  side  of  the  wafer  was  limited  in  extent  because  of  the  very  small  spot  size  of  the 
input  objective. 

The  SUT  was  then  aligned  about  the  x  axis  (Fig.  5-2)  by  tilting  the  sample  until  the 
film-coated  edge  of  the  wafer  could  be  seen  on  both  monitors  and  a  "sunrise"  pattern  was 
seen  on  the  output  side  of  the  wafer.  The  spot  and  rings  in  the  "sunrise"  pattern  of  Fig.  6-1 
are  not  an  image  of  the  spot  on  the  input  side.  These  features  result  from  the  Fresnel 
diffraction  of  the  input  spot  by  the  film  edges.  This  pattern  above  a  waveguide  indicates 
that  the  beam  and  sample  are  aligned  such  that  light  will  be  injected  into  that  waveguide. 

The  sample  could  be  examined  by  translation  along  the  x  axis  until  waveguide  features  or 
waveguiding  through  the  film  were  observed.  The  SUT  was  aligned  in  the  x-z  plane  about 
the  y  axis  by  maximizing  the  amount  of  light  coupled  through  a  waveguide.  Typically  the 
alignment  process  was  iterative  in  nature,  tilting  the  sample  about  both  the  x  and  y  axis  until 
the  SUT  was  properly  aligned. 


Fig.  6-1.  The  "sunrise"  pattern  seen  on  the  output  monitor  from  the  diffraction  of  the  input 
light.  This  pattern  indicates  that  the  SUT  is  aligned  with  the  input  beam  and  that  light  will 
be  injected  into  the  guide  directly  below  the  spot. 
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After  the  sample  was  aligned,  the  PBS  on  the  output  side  could  be  replaced  with  a 
mirror  to  maximize  the  amount  of  light  collected  by  the  output  camera.  By  using  a  mirror 
in  this  position,  lower  input  beam  powers  could  be  used  while  still  producing  an 
observable  image  on  the  output  side.  With  high  input  powers,  the  PBS  could  be  put  in 
place  again,  protecting  the  output  camera  from  damage. 

Input  beam  power  was  varied  by  using  the  two  neutral  density  wheels.  Data 
collection  was  conducted  with  the  laser  in  the  mode-locked  configuration.  This  was  done 
principally  to  protect  the  polymer  film  from  excessive  heating  while  in  CW  mode.  The 
focusing  power  of  the  input  lens  made  it  possible  to  exceed  the  optical  damage  threshold  of 
the  polymer  film  with  the  laser  running  in  the  CW  configuration.  While  the  peak  power  of 
the  mode-locked  laser  was  much  greater  than  the  CW  power,  the  sample  was  only  exposed 
to  the  power  during  the  "on"  portion  of  the  pulse  (approximately  2  ps)  and  was  not 
illuminated  for  the  subsequent  13.3  ns. 

Images  from  the  cameras  were  captured  using  frame-grabbing  software  on  a  personal 
computer.  The  intensity  profile  of  a  line  drawn  across  the  image  could  also  be  obtained 
using  this  software. 


6-3 


VII.  Experimental  Results  and  Analysis 

In  this  chapter,  the  results  of  the  waveguide  optical  characterization  experiments  are 
presented.  Only  the  strip-loaded  waveguides  were  optically  characterized.  Guiding  of  the 
input  beam  was  observed  through  the  polyphenylene  film  in  an  area  of  the  wafer  that  did 
not  contain  waveguide  structures.  Guiding  through  waveguides  and  coupling  of  optical 
energy  from  one  waveguide  to  another  was  also  observed  in  the  sample  under  test  (SUT). 

In  all  experimental  observations,  the  laser  was  mode-locked  at  X0  =  0.9  pm,  with  the 
pulse  width  ~  2  ps,  and  the  electric  field  of  the  light  polarized  in  the  x  direction.  The  SUT 
had  a  length  in  the  z  direction  of  L  =  3  mm.  Because  of  problems  in  obtaining  a  consistent 
cleave  quality,  only  this  length  of  sample  was  successfully  characterized.  Based  on 
modeling  and  estimations  of  the  losses  in  the  optical  beam  train,  the  input  coupling 
efficiency  was  estimated  to  be  30%.  The  actual  input  efficiency  is  highly  dependent  on  the 
interface  cleave  quality  and  is  thus  difficult  to  quantify. 

In  all  pictures  of  the  output  side  of  the  sample  the  contrast  has  been  altered  slightly  to 
enhance  the  reproduction  quality  of  the  figure.  All  intensity  measurements  were  made  prior 
to  the  contrast  adjustments.  Intensity  profile  plots  in  the  x  direction  have  been  smoothed 
with  a  three-point  moving  average.  The  profile  plots  are  from  a  single  line  scan  of  one 
video  frame.  Because  of  the  difficulties  in  analyzing  the  general  case  of  a  multi-channel 
directional  coupler  composed  of  non-identical  waveguides,  the  results  presented  in  this 
chapter  are  qualitative  in  nature. 

Guiding  of  Light  Through  a  Polyphenylene  Slab  Waveguide 

The  guiding  of  light  through  polyphenylene  was  observed  in  a  slab  waveguide 
structure  (see  Fig.  7-1).  The  slab  waveguide  was  actually  a  region  of  the  strip-loaded 
waveguide  wafer  that  was  not  near  any  etched  waveguide  features.  The  length  of  the 
sample  was  L-  3  mm.  The  light  was  confined  in  the  y  direction  by  the  cover  and  substrate 
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and  allowed  to  diverge  in  the  x  direction.  Assuming  an  input  beam  with  a  Gaussian 
profile,  a  waist  of  2  pm,  and  X0  =  0.9  |im,  the  beam  would  diverge  to  a  waist  of  0.48  mm 

in  3  mm  of  travel  through  the  polyphenylene  guide  layer  (33:67).  Thus,  the  guided  input 
beam  would  spread  to  fill  the  field  of  view  in  the  x  direction  and  appear  as  a  line  on  the 
output  side  (see  Fig.  7-1).  The  changes  in  observed  intensity  along  the  x  direction  are 
probably  due  to  variation  in  the  output  edge  cleave  quality  for  that  direction.  These 
intensity  variations  could  also  be  due  to  defects  in  the  polymer  films  or  foreign  object 
contamination  of  the  output  side  edge.  Note  that  since  the  spot  size  of  the  output  objective 
was  =  3  pm,  the  guide  layer  (tg  ~  0.67  pm)  of  the  SUT  was  not  actually  imaged.  The 
image  on  the  output  side  of  the  experimental  apparatus  is  a  convolution  of  the  spot  size  of 
the  output  objective  and  the  profile  of  the  light  coming  from  the  guiding  layer.  A  measured 
intensity  profile  in  the  x  direction  of  the  guided  light  is  shown  in  Fig.  7-2a.  The  guide  is 
single-mode  in  the  y  direction. 


nn=  1  (air)  0.6  pm 


M  Polythiophene  EH  Polyphenylene  U  Silicon  dioxide 
(cover)  (guide)  (substrate) 


0.67  pm 
2.2  pm 


(a) 


(b) 

Fig.  7-1.  Guiding  of  light  through  the  polyphenylene  film.  In  (a),  a  diagram  of  the  slab 
waveguide  structure  is  shown.  In  (b),  the  light  band  near  the  center  of  the  picture  is  light 
(X0  =  0.9  pm)  being  guided  through  3  mm  of  polyphenylene  film. 
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Fig.  7-2.  Measured  (a)  and  modeled  (b)  intensity  profiles  for  a  polyphenylene  slab 
waveguide.  The  large  intensity  spike  in  (a)  near  y  =  -0.5  is  due  to  the  light  above  the  cover 
(see  Fig.  7-lb).  In  both  profiles  the  left  hand  part  of  the  profile  corresponds  to  the  cover 
region  while  the  right  hand  part  corresponds  to  the  substrate  layer.  In  (b)  the  y  =  0  point 
corresponds  to  the  center  of  the  modeled  guide  while  in  (a)  the  y  =  0  point  is  arbitrarily 
placed. 


In  agreement  with  modeling,  the  mode  profile  in  the  y  direction  is  single-mode  and 
asymmetric.  The  intensity  decays  more  rapidly  on  the  cover  side  of  the  guide  than  on  the 
substrate  side.  If  the  cover  layer  were  thick  enough  to  be  considered  semi-infinite  in  the  y 
direction,  the  opposite  would  be  true  since  the  index  difference  at  the  cover/guide  interface 
is  much  smaller  than  the  index  difference  at  the  substrate/guide  interface  (see  Figs.  7- la,  3- 
7a).  However,  in  the  slab  waveguide  region  of  the  sample,  the  cover  layer  is  only  0.6  |im 
thick  -  thin  enough  that  the  electric  field  profile  is  influenced  by  the  cover/air  interface.  The 
large  difference  in  index  between  the  cover  and  the  surrounding  air  causes  the  field  to  decay 
rapidly  on  the  cover  side  of  the  guide.  This  geometry  was  modeled  by  using  the  effective 
index  of  the  guide  in  a  slab  waveguide  model  with  air  as  the  cover  layer  (Fig  7-2b).  Note 
that  the  shapes  of  the  mode  profiles  appear  to  be  in  good  agreement. 
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Coupling  in  a  Multi-channel  Directional  Coupler 


In  a  multi-channel  directional  coupler,  such  as  the  groups  of  strip-loaded  guides 
fabricated  as  part  of  this  research,  the  light  traveling  in  a  waveguide  is  coupled  to  adjacent 
guides  on  both  sides.  Thus,  a  branching  or  cascading  of  the  light  through  the  guides  is  to 
be  expected  as  the  light  travels  in  the  z  direction  (see  Fig.  7-3). 


z  =  0 


z  =  L 


Fig.  7-3.  Flow  of  injected  power  in  a  multi-channel  directional  coupler.  The  individual 
guides  are  not  identical  but  are  coupled  to  both  adjacent  guides.  The  input  power,  initially 
launched  into  a  single  guide,  branches  out  to  the  other  guides  in  a  cascade  manner. 


Somekh  (26)  has  analyzed  a  multichannel  directional  coupler  for  the  special  case  of  an 
infinite  number  of  parallel,  identical  guides,  each  weakly-coupled  to  only  the  adjacent 
guides  (26:46-52).  Solving  the  coupled-mode  equation  for  this  special  case,  assuming  that 
light  is  injected  into  only  one  guide  at  z  =  0,  gives  the  power  flow  in  the  vth  guide  as 
(26:48) 


Pv(z)  =  7^(2xz)exp[-az] 


(7-1) 
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where  Py(z)  is  the  power  flow  in  the  vth  guide  (v  =  0  is  the  guide  the  light  is  initially 
launched  into),  Jy  is  the  Bessel  function  of  the  vth  order  (v  =  0,  1,  2  ...),  Kris  the  coupling 
coefficient  between  adjacent  guides,  z  is  the  direction  of  propagation,  and  a  is  the  loss 
coefficient  for  the  guides. 

Thus,  given  the  relative  output  intensities  of  the  guides  for  a  given  distance  z  =  L,  the 
coupling  coefficient,  k,  can  be  determined.  Unfortunately,  the  guides  fabricated  in  this 
research  did  not  lend  themselves  to  this  analysis.  The  analysis  assumes  an  infinite  number 
of  identical  guides  but  there  are  only  five  guides  of  a  given  size  within  a  group. 
Furthermore,  the  guides  are  not  identical,  both  by  design  and  because  of  the  rough 
sidewalls  of  the  guides,  and  thus  the  coupling  coefficient  between  the  guides  is  not  a 
constant.  Because  of  the  rough  sidewalls  and  differing  guide  dimensions,  there  is  a  "built- 
in"  phase  mismatch  between  adjacent  guides  that  varies  in  an  essentially  random  manner 
along  the  direction  of  propagation.  The  nature  of  the  waveguides  fabricated  during  this 
research  limits  the  analysis  to  qualitative  observations. 

Lateral  Confinement  of  Light  and  Coupling  Between  Guides 

The  confinement  of  light  in  the  lateral  (x)  direction,  indicating  the  guiding  of  light 
through  the  fabricated  waveguides,  was  observed  in  the  strip-loaded  waveguides. 

Coupling  of  light  into  the  s  =  2  pm  group  of  guides  on  a  L  =  3  mm  sample  was  successful. 
Furthermore,  coupling  between  adjacent  guides  was  observed. 

As  shown  in  Fig.  7-4,  the  guides  of  the  fabricated  multi-channel  directional  couplers 
are  numbered  starting  with  the  leftmost  guide.  Guides  G1  -  G5  have  W  =  10  pm,  G6  - 
G10  have  W  =  5  pm,  and  G1 1  -  G15  are  the  W  =  3  pm  guides. 

In  Fig.  7-4,  light  was  injected  into  guide  3  (G3,  W  =  10  pm)  on  the  input  side  of  the 
sample.  Recall  that  the  spot  size  of  the  input  objective  was  2  pm,  much  smaller  than  the 
lateral  dimension  of  the  W  =  10  pm  guides.  The  light  is  coupled  into  the  adjacent  W  =  10 
pm  guides  and,  to  a  much  lesser  extent,  into  the  IT  =  5  pm  guides.  In  Fig.  7-4,  the 
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average  power  at  the  power  meter  location  (Fig.  5-2)  was  0.1  W,  indicating  a  peak  power 
of  665  W  at  this  point  in  the  optics  train.  Using  30%  as  the  estimated  input  coupling 
efficiency,  the  peak  injected  power  is  ~  200  W. 


Polythiophene 

(cover) 


Polyphenylene 

(guide) 


Silicon  dioxide 
(substrate) 


Fig.  7-4.  Observed  confinement  of  light  in  the  lateral  direction  and  coupling  between 
adjacent  strip-loaded  guides.  The  nominal  separation  between  guides  is  s  =  2  pm.  The 
experimentally  observed  output  side  of  the  sample  is  shown  in  (a).  A  diagram  of  the  guide 
structure,  with  arrows  pointing  to  the  corresponding  guides,  is  shown  for  reference.  Light 
is  injected  into  guide  3  (G3)  on  the  input  side  of  the  sample.  As  shown  in  (b),  coupling 
occurs  to  the  other  guides. 
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It  should  be  noted  that  while  the  profiles  presented  in  this  chapter  represent  only  one 
line  scan  in  the  *  direction,  the  features  of  a  given  profile  persisted  from  frame  to  frame  and 
did  not  fluctuate.  The  intensity  peaks  and  dips  of  the  profiles  are  not  "noise"  but  rather  are 
indicative  of  the  spatial  distribution  of  the  light  within  the  waveguides.  Some  of  these 
features,  mainly  the  smaller  ones,  may  be  due  to  scattering  off  of  defects  in  the  film  or 
interface  cleave.  The  larger  features,  such  as  those  seen  in  G2  and  G3  of  Fig.  7-4,  are 
most  likely  evidence  of  the  modes  propagating  in  the  waveguides. 

In  Fig.  7-4,  the  W  =  10  pm  guides  seem  to  be  capable  of  supporting  the  first  odd 
mode  in  the  lateral  direction,  as  predicted  by  modeling  in  Chapter  ID.  The  fundamental 
mode  is  excited  in  the  guide  by  the  focused  input  beam;  conversion  to  the  first  odd  mode 
occurs  from  scattering  in  the  guide  due  to  material  defects  or  the  roughness  of  the  defining 
sidewalls.  Modeling  of  the  electric  field  profiles  in  the  guide  region  indicate  that  the  field 
near  the  rough  side  walls  is  small.  The  main  influence  of  the  rough  side  walls  is  to 
effectively  change  the  width  of  the  guide  along  the  direction  of  propagation.  These  changes 
in  guide  width  change  the  mode  profiles  along  the  length  of  the  guide,  leading  to  a  loss  of 
energy  as  the  light  propagates  down  the  guide.  The  effect  is  similar  to  sending  an  RF 
signal  through  a  length  of  connected  coaxial  cables,  each  cable  having  a  slightly  different 
characteristic  impedance.  Reflections  due  to  impedance  mismatch  result  in  losses  to  the  RF 
signal  as  it  propagates  along  the  coaxial  cables. 

Confinement  of  the  light  in  the  lateral  direction  and  coupling  between  guides  was  also 
observed  in  the  SUT  by  injecting  light  directly  into  the  W  =  5  pm  guides  of  the  s  =  2  pm 
group.  In  Fig.  7-5,  the  average  power  measured  at  the  power  meter  insertion  point  was  50 
mW,  indicating  a  peak  power  of  335  W  at  this  point  in  the  optics  train  and  a  peak  injected 
power  of  ~  100  W. 
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Fig.  7-5.  Observed  guiding  and  coupling  of  light  injected  into  a  W  =  5  pm  strip-loaded 
waveguide.  The  nominal  separation  between  guides  is  s  =  2  pm.  The  view  is  of  the  output 
side  of  the  sample.  The  experimentally  observed  output  side  of  the  sample  is  shown  in  (a). 
A  diagram  of  the  guide  structure,  with  arrows  pointing  to  the  corresponding  guides,  is 
shown  for  reference.  As  seen  in  (b),  light  injected  into  G9  on  the  input  side  of  the  sample 
couples  into  the  adjacent  guides  in  a  cascade  or  branching  manner. 


No  light  was  observed  out  of  the  W  =  3  jam  waveguides.  This  is  probably  due  to  the 
rough  sidewalls  of  the  guides.  Because  the  nominal  width  of  these  guides  is  smaller,  the 
roughness  of  the  sidewalls  would  have  a  more  drastic  effect  on  the  effective  width  of  the 
guide,  resulting  in  greater  losses.  Any  light  coupled  into  a  W  =  3  pm  guide  is  probably 
lost  before  reaching  the  output. 

In  agreement  with  modeling,  the  W  =  5  pm  guides  G8  and  G9  appear  to  be  single¬ 
mode  in  the  x  direction.  The  first  odd  mode  appears  to  be  excited  in  G10  and  G7.  Flaws 
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along  the  length  of  the  guides  could  change  the  effective  width  of  the  guides  and  allow  the 
first  odd  mode  to  propagate  for  some  distance  along  G10  and  G7.  These  intensity  profiles 
could  also  be  due  to  contamination  of  the  output  coupler  edge  at  these  guide  locations. 

The  relatively  large  amount  of  coupling  from  G6,  a  W  =  5  fim  guide,  into  G5,  a  W  = 
10  pm,  is  to  be  expected.  Modeling  indicates  that  the  first  odd  mode  is  poorly  confined  in 
the  W  =  10  |im  guides.  Thus  the  first  odd  mode  in  G5  would  be  well  coupled  to  the 
fundamental  mode  of  G6.  In  Fig.  7-5a  the  light  is  injected  into  G9  rather  than  G8,  the 
center  guide  of  the  W  =  5  (im  group,  in  an  attempt  to  better  isolate  the  W  =  5  p.m  guides 
from  the  W  =  10  |im  guides. 

Another  example  of  coupling  between  the  waveguides  is  shown  in  Fig.  7-6.  Here, 
light  is  injected  on  the  input  side  into  G8.  The  average  power  measured  at  the  power  meter 
insertion  point  was  50  mW,  implying  a  peak  injected  power  of  «  100  W.  The  intensity 
profiles  of  guides  G6  -  G9  in  Fig.  7-6  look  similar  to  the  corresponding  profiles  in  Fig.  7- 
5.  In  both  figures  the  guide  into  which  light  was  initially  launched  has  the  highest  relative 
intensity  on  the  output  side  of  the  sample.  This  is  probably  due  to  the  built-in  phase 
mismatch  between  the  guides,  which  decreases  the  coupling  between  guides.  It  is  also 
possible,  given  the  asymmetric  and  non-identical  nature  of  the  waveguides,  that  a  majority 
of  the  power  has  returned  to  G8  at  this  particular  sample  length.  The  problems  associated 
with  waveguide  fabrication  and  cleaving  make  a  complete  analysis  of  this  data  difficult. 

Note  that  in  both  G4  and  G5  of  Fig.  7-6a,  the  intensity  profile  is  not  centered  on  the 
waveguide,  but  rather  is  shifted  spatially  to  one  side.  This  is  probably  the  result  of  both  the 
fundamental  and  first  odd  modes  propagating  in  the  waveguide.  If  more  than  one  mode  is 
propagating  in  a  waveguide,  the  electric  field  in  the  guide  will  be  the  sum  of  the  electric 
fields  of  the  modes.  Within  the  guide  there  were  will  be  regions  where  the  electric  fields 
add  in  such  a  way  as  to  cancel  each  other.  Likewise,  there  will  be  regions  where  the 
electric  fields  of  the  modes  add  and  reinforce  each  other.  Thus,  it  is  possible  that  the 
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fundamental  and  first  odd  mode  are  propagating  in  G4  and  G5,  resulting  in  an  overall 
intensity  pattern  as  shown.  Note  that  both  of  these  guides  seem  to  have  good  output 
cleaves,  as  evidenced  in  Fig.  7-5.  This  suggests  that  the  intensity  patterns  of  G4  and  G5  in 
Fig.  7-6a  are  not  due  to  interface  contamination  or  poor  quality  of  the  output  cleave. 
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Fig.  7-6.  Coupling  of  light  between  adjacent  strip-loaded  waveguides.  The  nominal 
separation  between  guides  is  s  =  2  pm.  The  experimentally  observed  output  side  of  the 
sample  is  shown  in  (a),  along  with  a  diagram  of  the  guide  structure.  As  seen  in  (b),  light 
injected  into  G8  on  the  input  side  of  the  sample  couples  into  the  adjacent  guides. 
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Nonlinear  Effects 


As  the  power  injected  into  a  multi-channel  NLDC  is  increased  from  a  low  level  to  a 
level  large  enough  to  induce  nonlinear  effects,  some  variation  of  the  coupling  coefficient, 
k,  and  thus  the  output  intensities  would  be  expected.  To  investigate  the  nonlinear 
characteristics  of  the  fabricated  structures,  the  peak  power  injected  into  the  SUT  was  varied 
and  the  output  monitored.  Using  30%  as  the  input  coupling  efficiency,  the  estimated 
injected  power  varied  from  ~  41  W  to  «  319  W. 

The  critical  power  calculated  from  eqn.  (3-8)  does  not  accurately  characterize  the 
nonlinear  response  of  the  fabricated  multi-channel  NLDC  because  it  was  developed  by 
assuming  a  pair  of  identical,  isolated  nonlinear  waveguides  that  had  an  interaction  length  of 
L0.  The  phase-mismatch  of  the  guides  makes  the  modeling  much  more  difficult  (13:1581). 
At  best,  the  critical  powers  for  the  strip-loaded  guides  found  in  Table  3-1  can  be  regarded 
as  only  an  order  of  magnitude  estimate  of  the  power  needed  to  induce  nonlinear  effects  in 
the  fabricated  structure.  For  reference,  the  critical  power  for  a  pair  of  identical,  isolated 
strip-loaded  guides  with  W  =  5  (im  and  s  =  2  pm  was  estimated  to  be  230  W. 

In  Fig.  7-7,  the  output  profiles  for  two  different  injection  powers  is  shown.  It  should 
be  noted  that  the  profiles  represent  only  one  scan  of  the  output  intensities  and  are  thus 
essentially  qualitative  in  nature.  Note  that  the  basic  features  of  the  two  profiles  in  Fig.  7-7 
are  the  same.  This  fact  is  further  evidence  that  the  profile  features  are  representative  of  the 
spatial  distribution  of  light  within  the  guides  and  not  "noise."  In  both  profiles  the  power 
was  initially  launched  into  G9,  just  as  in  Fig.  7-5. 

As  can  be  seen  in  Fig.  7-7,  there  is  some  shift  in  the  intensity  profiles  for  the  two 
different  input  powers.  In  the  case  of  the  high  injection  power,  the  highest  peak  intensity 
was  observed  out  of  G8,  even  though  the  power  was  initially  launched  into  G9.  This  shift 
in  peak  output  intensity  from  G9  to  G8  could  indicate  that  the  coupling  between  G9  and  the 
adjacent  guides  changed,  as  would  be  expected  for  an  NLDC. 
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Fig.  7-7.  Output  intensities  for  "low"  and  "high"  peak  injected  powers.  In  both  profiles, 
all  of  the  light  was  initially  injected  into  G9.  The  shifts  in  the  output  intensities  could  be 
indicative  of  a  nonlinear  response  but  the  results  are  not  conclusive. 


Arguably,  there  may  be  significant  shifts  in  output  intensity  for  G4,  G5,  and  G6, 
possibly  indicating  a  change  in  coupling  due  to  higher  injection  power.  However,  without 
a  more  quantitative  analysis,  which  is  made  difficult  by  the  multi-channel,  non-identical 
nature  of  the  fabricated  guides,  no  definitive  conclusions  can  be  reached  concerning  the 
nonlinear  response  of  the  SUT.  The  observation  of  nonlinear  effects  would  be  more 
straight-forward  in  devices  consisting  of  a  single  pair  of  isolated,  identical,  single-mode 
waveguides  since  the  influence  of  multiple,  nonidentical  adjacent  guides  would  not  be  a 
factor  in  the  analysis.  A  much  more  quantitative  investigation  of  nonlinear  effects  on 
coupling  would  be  possible  with  such  devices. 

Waveguide  Losses 

Measurements  of  the  waveguide  loss  were  not  made  due  to  the  difficulty  in  obtaining 
consistent  cleave  quality  as  the  sample  length  was  decreased.  The  analysis  of  loss 
measurements  would  have  been  further  complicated  by  the  presence  of  multiple  coupled 
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guides.  Measurements  could  have  been  made  on  the  channel  waveguide  structures  since 
they  were  essentially  uncoupled  but  these  measurements  would  not  have  been  applicable  to 
the  analysis  of  the  strip-loaded  waveguides.  Also,  the  poor  polymer  quality  observed  in 
the  channel  waveguide  fabrication  (Fig.  4-11)  indicated  that  the  film  may  not  have  had 
good  optical  qualities.  The  scattering  losses  from  this  inhomogenous  film  would  have 
probably  dominated  over  the  scattering  from  the  substrate  and  air  interfaces. 
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VIII.  Summary  and  Recommendations 

In  summary,  this  thesis  research  has  resulted  in  the  successful  design  and  fabrication 
of  waveguiding  structures  in  polyphenylene,  a  recently  developed  nonlinear  organic 
polymer.  Waveguiding,  through  both  a  slab  waveguide  and  strip-loaded  guides,  was 
observed  for  the  first  time  in  this  polymer.  The  coupling  of  optical  energy  from  one  guide 
to  the  adjacent  guides  in  a  multi-channel  directional  coupler  was  also  observed.  No 
definitive  conclusions  concerning  nonlinear  effects  are  possible  due  to  the  non-identical, 
multi-channel  nature  of  the  fabricated  waveguides.  Several  design  and  modeling  tools  have 
been  developed  that  should  support  additional  research  in  this  area. 

The  film  growth  results  for  the  strip-loaded  guides  indicate  that  it  is  possible  to  grow 
both  polyphenylene  and  polythiophene  films  that  are  dense,  uniform,  mechanically  and 
hydroscopically  stable,  and  suitable  for  fabrication  into  IOCs.  However,  the  film  deposited 
for  the  channel  waveguides  was  not  of  the  same  quality.  Reproducibility  of  film  quality 
should  be  quantified  with  further  study.  Better  attention  to  substrate  roughness  and  the 
history  of  the  distilled  benzene  or  dichlorothiophene  should  help  improve  film  quality. 

Another  difficulty  encountered  in  the  fabrication  of  waveguides  involved  the  growth 
rate  of  the  thin  polymer  films.  Slow  film  deposition  rates  hinder  the  fabrication  of 
waveguides.  In-situ  monitoring  of  film  thickness  would  greatly  enhance  the  ability  to  grow 
films  to  a  specified  thickness. 

One  way  to  avoid  the  CVD  growth  of  another  polymer  as  a  cover  layer  (such  as  the 
use  of  polythiophene  in  the  strip-loaded  guides)  would  be  to  spin-on  a  layer  of  photoresist 
direcdy  on  top  of  the  polyphenylene  and  use  it  as  the  cover  layer.  Most  common 
photoresists  have  refractive  indices  of  ~  1.6  and  are  essentially  transparent  to  near  infrared 
radiation,  making  them  suitable  materials  for  such  an  application.  The  photoresist  could 
simply  take  the  place  of  the  polythiophene  layer  in  Fig.  4-2,  be  passivated  with  silicon 
nitride,  and  etched  with  oxygen  in  a  REE  chamber.  Alternatively,  the  photoresist  on  top  of 
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the  polyphenylene  could  be  exposed  with  a  lithographic  mask  and  be  developed,  leaving 
the  polyphenylene  exposed.  The  large  step  in  index  at  the  polyphenylene  /air  interface 
would,  however,  tend  to  increase  the  confinement  of  light  in  the  lateral  direction, 
decreasing  the  coupling  between  guides  and  the  maximum  allowable  guide  width  for 
single-mode  operation. 

Another  concern  with  using  photoresist  as  a  cover  layer  is  its  stability  when  exposed 
to  UV.  Unless  the  photoresist  was  passivated  with  silicon  nitride  or  a  similar  material  after 
waveguide  fabrication,  the  photoresist  unexposed  during  fabrication  could  become  exposed 
(thus  changing  its  optical  and  mechanical  properties)  during  subsequent  handling  and 
testing.  The  potential  benefit  of  using  photoresist  as  the  cover  layer  -  avoiding  the  slow 
CVD  growth  of  another  polymer  on  top  of  the  polyphenylene  -  justifies  an  attempt  at  this 
fabrication  method.  The  use  of  a  photoresist  cover  layer  should  be  considered  during 
further  research. 

The  results  of  this  thesis,  the  successful  fabrication  of  guides  and  the  observation  of 
both  guiding  and  coupling  of  light,  justify  further  investigation  of  polyphenylene  as  a 
material  for  IOCs.  Perhaps  the  most  important  next  step  in  this  research  is  the  fabrication 
of  a  photolithographic  mask  more  suitable  for  NLDC  research.  The  mask  should  contain 
pairs  of  single-mode  waveguides,  coupled  to  each  other  but  far  enough  away  from  other 
waveguides  to  be  effectively  isolated.  NLDCs  with  a  variety  of  guide  dimensions, 
separations,  and  coupling  lengths  could  be  incorporated  in  the  mask,  allowing  the 
researcher  to  more  fully  characterize  the  coupling  between  guides  and  experimentally 
confirm  the  nonlinear  optical  properties  of  polyphenylene.  Isolated  waveguides,  for 
making  loss  measurements,  could  also  be  incorporated  in  the  mask  design.  An  example  of 
such  a  mask  is  shown  in  Fig.  8-1  . 
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Fig.  8-1.  Proposed  lithographic  mask  for  further  NLDC  research.  The  registration  marks  serve  to  identify  the  proper 
cleave  locations  and  the  bar  channels  of  the  devices.  The  separation  between  couplers  should  be  large  (~  20  pm)  so 
that  adjacent  couplers  are  isolated.  All  the  couplers  in  a  column  could  be  identical  while  the  couplers  along  a  row 
could  have  different  guide  separations  and  interaction  lengths. 


Appendix  A:  Waveguide  Dimensions  for  Single-mode  Operation 

This  document  calculates  the  min  and  max  waveguide  dimensions  for  a  single-mode 
waveguide.  The  indicies  of  the  materials  and  the  free-space  wavelength  of  operation  must 
be  entered  into  the  document.  See  Fig.  2-2  for  geometry. 
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Appendix  B:  Slab  Waveguide  Electric  Field  Profiles 


This  Mathcad  document  calculates  and  plots  the  profile  of  the  mode  in  a  slab  waveguide. 
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mode  ■  =  0  Mode  to  calculate  profle.  Note  that  m=0  is  the  fundamental  mode.  If  a  mode  is 
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Appendix  C: 

Effective  Index  Method  for  Strip-Loaded  Waveguides 

This  document  uses  the  effective  index  method  to  calculate  the  effective  index  of  Regions  I  and  II  for 
the  strip-loaded  waveguide  structure  shown  below. 
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Calculation  of  the  Drona&ation  constant  in  the  z  direction 


Trancendental  equations  to  be  solved  for  beta.  Note:  absolute  value  of  argument  because  of  1st  quadrant! 
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Finding  beta  for  region  I  (outside  loading  strip). 
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Appendix  D: 

Effective  Index  Method  for  Rib  Waveguides 


This  document  uses  the  effective  index  method  to  calculate  the  effective  index  of  Regions  I  and  II  for 
the  rib  waveguide  structure  shown  below. 
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Calculation  of  the  propagation  constant  in  the  z  direction 
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Trancendental  equations  to  be  solved  for  beta.  Note:  absolute  value  of  argument  because  of  1st  quadrant! 
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Finding  beta  for  region  II  (under  the  loading  strip). 
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NOTE:  The  initial  guess  is  critical  to  the  success  of  this  routine.  Use 
the  above  graph  to  pick  a  beta  that  is  very  close  to  the  value  at  intersection. 
If  Mathcad  has  a  problem  finding  a  solution  that  obviously  exists,  try 
increasing  the  tolerance  of  the  Diff(Beta)  function  in  the  Given  block. 


P>  n  0  k  o 

RegionIIDiff(P)<  0.0001 

Regionlip  m  !=  fmd(  P) 

RegionIIpm  =  1.115*  107’m  1 
Regionlip  m 

RegionllN  eg-  != -  RegionllN  eff  =  1.632 

k  Q 


p  :=  1.114-  107  m_1 
given 

P<  n  g'k  o 
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Finding  beta  for  region  I  (outside  loading  strip). 

P  ■  -  n  o'  k  o>  1.05' n  q  k  q. .  n  g  k  ^ 


RegionIDiff(  (1)  :=  |RegionILHS(  P)  -  RegionIRHS(  P)  | 

NOTE:  The  initial  guess  is  critical  to  the  success  of  this  routine.  Use 
the  above  graph  to  pick  a  beta  that  is  very  close  to  the  value  at  intersection. 
If  Mathcad  has  a  problem  finding  a  solution  that  obviously  exists,  try 
increasing  the  tolerance  of  the  Diff(Beta)  function  in  the  Given  block. 


P>  n  o  k  o 


P  -=  1.09-  107  m~ 1 
given 

P<  n  g  k  o 


RegionIDifff  p><  0.0001 
Regionip  m  !=  find(  p ) 

Regionip  m  =  1.101*  107,m_1 

mode  y  =  0 

Mode  number  in  the  x  direction 
for  which  these  numbers  were 
calculated.  They  ARE  mode 
dependent!! 

N  eq  '=  (n  g2  —  ReSionIIN  eff2  +  Re§ionIN  eff2) 


RegionINeff  - 


Regionip 


m 


RegionIN  =  1.612 

RegionllN  eff  =  1.632 
1 


Naj=  1.691 


Note:  This  is  the  value  that  should  be  substituted  for  ng  and  05 
in  Maracatilli's  method. 
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Appendix  E:  Marcatili's  Method  for  Embedded  Guides 

This  MathCAD  document  calculates  and  plots  the  profile  of  the  transverse  modes  in  a  general 
waveguide  structure  like  the  one  shown  below. 


Constants 


pm 

:=  10_6m 

nm 

:=  io"9m 

c  :=  2.998- 108*  — 

sec 

kW 

:=  103*watt 

GW  : 

=  109*  watt 

“3 

gram  ■-  10  -kg 

Material  and  Structure  Parameters 


n  J  :=  1.79 
n  2  :=  1-71 
n  3  :=  1.789 
n4  :=  1.5 
n  5  :=  1.789 


W  :=  5- pm 

X  0  :=  0.92- pm 

i  —2'n  —i 

ko  -7—  ®  -koc 

A,  o 

t  g  :=  0.67- pm 

Area  -=  t  •  W 

O 

s  :=  2- pm  guide  separation  in  coupler 

f  w 
aspect  • - 

lg 

aspect  =  7.463 

Note:  Method  works  best  for  aspect  ratios  greater  than  2, 

p  :  =  1  for  modes  in  x  direction 
q  :  =  1  for  modes  in  y  direction 

Mode  numbers.  Note  that  p  =  q  =1  is  the  fundamental  mode. 


2mks  :=  2.03- 10' 10 

<-  NOTE:  n2  MUST  be  in  units  of  cmA2/kW 

This  is  the  value  for  polyphenylene. 

n  2mks'n  1 

sec2 

pm  • 

2esu  '  g  4 

2esu  =  4.326*  10~ 1 1 

kill 

gram 

sec2 

’em - 

gram 

See  formula  A5.13b  from  (3),  pg  306-308 

These  are  tlie  ESU  units  used  by  Jensen, 
n  =  no+n2esu*IEIA2  so  n2esu  has  units  of  1/EA2  in 
ESU  units:  cmA2/statvoltA2  =  cm*secA2/gram 
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LHS  y(  k  y )  :=  k  y  t  g  Equations  to  solve  for  propagation  constants,  from  Marcatili. 

RHS  y(ky)  :=qjt-atan(  Jk  y-r|  2< k  y )  |  >  “  atan(  |k  yri  4(  k  y )  |  ) 
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Diffy(ky)  :=  |LHSy(ky)-RHS  y(ky)| 

k  :=  3.5- 106’  m  1  NOTE:  The  initial  guess  is  critical  to  the  success  of  this  routine.  Use 

the  above  graph  to  pick  a  beta  that  is  very  close  to  the  value  at  intersection. 

given 

ky>0m  1 

k  yCk  o 

Diffy(ky)<  0.0001 

k  y  :=  find(  k  y )  k  y  =  2.758*  106  *m  1  Propagation  constant  in  y  direction 


k  v  :=  0-  m”  \  20000  m"1..  0.5- kv 

a  y 


Diffx(kx)  :=  LHSx(kx)-RHSx(kx) 


k  4- 105-  m  1  NOTE:  The  initial  guess  is  critical  to  the  success  of  this  routine.  Use 

the  above  graph  to  pick  a  beta  that  is  very  close  to  the  value  at  intersection. 


given 

kx>0m  1 
kx<k0 

Diffx(kx)<  0.0001 

kx:=find(kx)  k  x  =  2.995*  105  *m” 1 


Propagation  constant  in  x  direction 
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Plot  profile  of  mode  in  v  direction 


R2  G1( y>Cy)  :=Cy  exp 


■n2(ky) 


R2Gj  =  region  2  for  guide  1 
RlyGl  =  region  1  in  y  for  guide  1 
(follows  this  pattern,  see  figure  aat 
beginning  of  document  for  region 
locations.) 


Riy  gi^  y>  c  y)  ,=  c  y- 


cos  k , 


— 1 — Vsi„ft  (y_ls)  1 

V  2/  J  yk y-ri  2(ky)y  L  y  v  2/ 


R4 


Gl(y’Cy>  :=C  y 


COS(  ky‘("tg))  — 


(  '  Vs, 

\k  y'Tl  2(ky  )/ 


sin(ky-(-tg)) 


I  t.1 

fy+-r 

•exp 

2 

^^ky) 

RlR2G1(y,Cy)  :=if|y>-^,R2G1(y,Cy),RlyG1(y,Cy)j 
yeprofile(  y,C  y)  :=it(y<-^,R4G1(y,  C  y).RlR2  G1(y,C  y)) 


Normalization  y(  C  y ) 


3t 


"3-t , 


yeprofde(y,C  yY 


dy 


See  Saleh,  pg.  251 


C  y  :=  960  m 


-0.5 


initial  guess  for  normalization  constant 


given 

r  ^  n  “0.5 
L  y>  Om 

Noimalization  y(  C  ) = 1 
Cy  :=fmd(Cy) 

C  y  =  960.383* m  0  5  Normalization  y(C  )  =  1 
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Plotting  the  electric  field  profile 

y  •=  "  1 .5- 1  g ,  “  1 .45- 1  g  •  •  2.5- 1  g  Electric  Field  Profile 


.  _n*7*  Note:  y  =  0  is  center  of  guide.  Positive  y  values  are 

guide  thickness  is.  t  g  -  U.6/  ftm  Upper  half  of  guide  and  cover.  Negative  y  values  are  lower 

half  of  guide  and  substrate. 

mode  number  in  y  direction  is:  q  =  1  Note:  q  =  1  is  fundamental  mode 


Plotting  profile  of  mode  in  x  direction 


R3  Gf(x,C  x)  :=C  x  exp 


M 

S3(kx) 


RlxG1(x,Cx)  :=Cx- 


cos  k 


R5  G1(x,C  x)  :=C  x-  cos(kx  (-W)) 


(kx^3(kx)) 


•sin(k  X-(“W) )  exp 


x  + 


1*5  (kx>j 


R1R3  G1(  x,  C  x )  :=  if  ^x>  ^ ,  R3  G1(  x,  C  x ) ,  Rlx  G1(  x,  C  x  )j 

/  W  \ 

xeprofilel(  x,  C  x)  •=  iflx<~— ,  R5  G^(  x,  C  x),  R1R3  Gj(  x,  C  X)1 
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to|  ^ 


1 2-W 

Normalization  x(  C  x )  xeprofilel(  x,  C  x  )2  dx  from  Saleh 

. -2-W 


C  x  •“  180' m  0-5  initial  guess  for  normalization  constant 
given 

Cx>0m"°-5 
Normalization  x(  C  ) = 1 

Cx  :=fmd(C  x) 

C  x  =  297.605* m  0  5  Normalization  X(C  x)  =  1 
Plotting  the  electric  field  profile 

x  :=-2W,-1.9W..  2-W  Electric  Field  Profile 


guide  width  is:  W  =  5  *pm  Note:  x  =  0  is  center  of  guide.  Positive  x  values  are  right 

half  of  guide  and  cover.  Negative  x  values  are  left  half  of 
guide  and  substrate,  (see  figure  at  beginning  of  document) 

mode  number  in  x  direction  is:  p  =  1  Note:  p  =  1  is  fundamental  mode. 
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Calculation  of  coupling  between  guides 


Field  profile  for  second  guide,  centered  on  x  =  s  +  W 


R3  G2(x)  :=Cx  exp 


(-?-—) 

^3(*x) 


R1xG2(x)  :=Cx- 


cos  k 


.(5_|_s_w)]_(_^).^kx.(x_w_s_v 


RlxR3  Q2^  x)  :=  if^x>  ^s  +  }W)’R3G2(x)’R1xG2(x 
xeprofile2(  x )  :=  Lf|x<  ^s  +  -i- wj  ,  R5  G2<  x ) ,  RlxR3  G2( 


•sin(k  -(-W))  -exp 


W 

x  + - s- W 

2 

^5(kx> 


x) 


xIprofile2(  x )  1=  xeprofile2(  x ) 
x  :=- 5- jxm, -4.9- (Am..  15  |xm 


Electric  Field  Profiles 


guide  width  is:  W  =  5  "Jim 
guide  separation  is:  s  =  2  •ji.m 
Note:  p  =  1  is  the  fundamental  mode. 


Note:  x  =  0  is  the  center  of  guide  1. 
mode  number  in  x  direction:  p  =  1 
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Closely  Coupled  Calculation  Marcatili,  pg.  2075,  2099 

N  eq  :=  n  3  k  z  :=  (k  o2’n  l2  “  k  x2  “  k  y2)  2  k  z  =  U91’  107-m_1 

kxeq:=[<koN<q>2-kz2-ky2]2 

|sk  xeq|  =  0.556  The  guides  are  weakly  coupled  if  this  quantity  is  »1 . 


2'k>-%3<kx)"seXPl 

f-  s  \ 

i  5  3 <k*V 

kzl 

(tg+sin(kytg)2tl2(ky))(w+2^3(kx))' 

(k  x  +  2) 

\  53(kx)2/ 

k  m  =  0.741  'mm 


■  w 

k  := 

k  2 

1  r  2  2\  ko 

2  '  1  n3jl“ 

2 

xeprofilel(  x,  C  x )•  xeprofile2(  x )  dx 

Using  Saleh’s  equation 

z 

_W 

< 

2 

k  =  1.185  ’mm  1  W  =  5  ‘pm  s  =  2  ‘pm 

p  =  1  mode  number  in  x  direction 


L  „  := — — —  L  „  =  1.326  'mm 
0  (2-k)  0 


Calculation  of  Critical  Power 


^  o' ®  l'c 

P  „T  :=  k-  Area - 

2’n2’n  2esu 


P  cj  =  229.51  ’watt 
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